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ABSTRACT 

We compare molecular gas traced by "'^^00(2-1) maps from the HERACLES survey, with tracers of 
the recent star formation rate (SFR) across 30 nearby disk galaxies. We demonstrate a first-order linear 
correspondence between Smoi and EgpR but also find important second-order systematic variations 
in the apparent molecular gas depletion time, t^"' = Snioi/^sFR- At the 1 kpc common resolution of 
HERACLES, CO emission correlates closely with many tracers of the recent SFR. Weighting each line 
of sight equally, using a fixed alphaco equivalent to the Milky Way value, our data yield a molecular 
gas depletion time, r"^"' = Emoi/SgFR « 2.2 Gyr with 0.3 dex la scatter, in very good agreement 
with recent literature data. We apply a forward-modeling approach to constrain the power-law index, 
A'', that relates the SFR surface density and the molecular gas surface density, "Ssfr oc S^^j. We find 
A'' = 1 ± 0.15 for our full data set with some scatter from galaxy to galaxy. This also agrees with 
recent work, but we caution that a power law treatment oversimplifies the topic given that we observe 
correlations between r^^"' and other local and global quantities. The strongest of these are a decreased 

T^p in low-mass, low-metallicity galaxies and a correlation of the kpc-scale r^""' with dust-to-gas ratio, 
D/G. These correlations can be explained by a CO-to-H2 conversion factor (aco) that depends on 
dust shielding, and thus D/G, in the theoretically expected way. This is not a unique interpretation, 
but external evidence of conversion factor variations makes this the most conservative explanation of 
the strongest observed t^'^ trends. After applying a D/G-dependent aco, some weak correlations 

between t^"^ and local conditions persist. In particular, we observe lower t^'^ and enhanced CO 
excitation associated with nuclear gas concentrations in a subset of our targets. These appear to 
reflect real enhancements in the rate of star formation per unit gas and although the distribution of 
Tdcp does not appear bimodal in galaxy centers, Tdop does appear multivalued at fixed Smoi, supporting 
the the idea of "disk" and "starburst" modes driven by other environmental parameters. 
Subject headings: galaxies: evolution — galaxies: ISM — radio lines: galaxies — stars: formation 



1. INTRODUCTION 

The relationship between gas and star formation 
in galaxies plays a key role in many areas of astro- 
physics. Its (non)evolution over cosmic time informs 
our understanding of gala xy evolu t ion at hig h rcdshift 
dPadd i ct al. 2010; Tacco ni et all 120101 : l^nzcl ct aD 
The small-scale efficiency of star formation is 
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a key input to galaxy simulations and scaling relations 
measured for whole galaxies provide important bench- 
marks for the output of these simulations. Measurements 
of gas and star formation at large scales give context 
for studies focu sing on parts of the Milky Way (e.g., 
iLada et "all 1201 0: Hcidc rman et a l. 2010) and the near - 
est galaxies fe.g.. iSchruba et al.ll2010l: iChen et al.ll2010D . 
Ultimately, a quantitative understanding of the gas-stars 
cycle is needed to understand galaxy evolution, with im- 
plications for the galaxy luminosity function, the galaxy 
color-magnitude diagram, the structure of stellar disks, 
and chemical enrichment among other key topics. 

Recent multiwavelength surveys make it possible to es- 
timate the surface densities of gas and recent star forma- 
tion in dozens of nearby galaxies. This has lead to several 
studies of the relationship between ga s and stars. Many 
of these focus on a sing le galaxy (e.g.. IHeyer et al.112004: 
Kennicutt et al.l l2007t iBlanc et all 120091: iVerlev et al.l 
2010': 'Rahman et al.' '2 011D or a sma ll sample (e.g., 
Wilson ct al. 2009; War ren et aIll2010D . Restricted by 
the availability of complete molecul ar gas maps, stud- 



ies of larg e sets of galaxie s (e.g 
Kennicutt ll998bl :_[Rownd fc Youm 



l2002i,Lerov et al 



lYoung et all 119961: 

19991: iMurgia et"all 

[2005uSaintonge et al.ll201lD mostlv use 



integrated measurements or a few low-resolution point- 
ings per galaxy. 
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From 2007-2010, the HERA CO- Line Extragalactic 
Survey (HERACLES, first maps in iLerov et alj [20091) 
used tlie IRAM 30-m telescopeQ to construct maps 
of CO emission from 48 nearby galaxies. The com- 
mon spatial resolution of the survey is ^ 1 kpc, suf- 
ficient to place many resolution elements across a typ- 
ical disk galaxy. Because the targets overlap surveys 
by Spitzer (mostly SINGS and LVL, iKennicutt et al.l 
2003a; Dale et al. 2009) and GALEX (mostly the NGS, 
Gil de Paz et al.l 120071) . a wide variety of multiwave- 
length data are available for most targets. In this paper, 
we take advantage of these data to compare tracers of 
molecular gas and recent star formation at 1 kpc resolu- 
tion across a large sample of 30 galaxies. 

This paper buil ds on wor k bv | Lerov et al.l (|2008l here- 
after LOS) and Big iel et al.l (|2008|, hereafter BOS). They 
combined the first HERACLES maps wi th data from 
The H I Nearby Galaxies Survey (THINGS IWalter et al.l 
[2008I) . SINGS, and the GALEX NGS data to com- 
pare H I, CO, and tracers of recent star formation in 
a sample of nearby galaxies. In the disks of large spi- 
ral galaxies, they found little or no dependence of the 
star formation rate per unit molecular gas on environ- 
ment. The fraction of interstellar gas in the molecular 
phase, on the other hand, varies strongly within and 
among galaxies, exhibiting correlations with interstel- 
lar pressure, stellar surface density, an d to tal gas sur- 
face density among other q uantities (Wong_&: Blitz 2002;; 
iBlitz fc Rosolowskvl[200l LOS). They advocated a sce- 
nario for star formation in disk galaxies in which star 
formation in isolated giant molecular clouds is a fairly 
universal process while the formation of these clouds out 
of the atomic gas reservoir depends sensitively on envi- 
ronmcnt (see also, Wong 2009: Ostri ker et aLl(2010D . 

The full HERACLES sample spans a wider range of 
masses, morphologies, metallicities, and star formation 
rates (SFRs) than t he sp irals studied in L0 8 and BOS. 
iSchruba et all (|20lH l and lBigiel eFall (|20lH ) used these 
to extend the fin dings of LOS and BOS. Using stacking 
techniques, Schru ba et al.l ([2011) demonstrated that cor- 
relations between star formation tracers and CO emission 
extend into the regime where atomic gas dominates the 
ISM, Ehi > Sjiioi. This provides the strongest evidence 
yet that star formation in disk galaxies can be separated 
into star formation from molecular gas and the balance 
between atomic and i nolecular gas, a hypothe sis that has 
a long history (e.g.. lYoung fc Scovilld [l99ll and refer- 
ences therein). iBigiel et al.l ()201lD demonstrated that a 
fixed ratio of CO to recent star formation rate remains a 
reasonable description of the ensemble of 30 galaxies. 

In t his p aper we expand on LOS, BOS, IBigiel et al.l 
(poll . and lSchrubaetall ([20ll and examine the gen- 
eral relationship between molecular gas and SFR in 
nearby disk galaxies. We divide the analysis into two 
main parts. In Section [3l we consider the scaling rela- 
tion linking SgpR to Smoi, the molecular analog to the 
"Kennicutt-Schmidt law" or "star formation law." We 
show the distribution of data in SspR-Snioi parameter 
space using different weightings (Section 13.11) and ex- 
amine how this distribution changes with different ap- 
proaches to physical parameter — varying the choice of 

12 IRAM is supported by CNRS/INSU (France), the MPG (Ger- 
many) and the ION (Spain). 



SFR tracer, the CO transition studied, the processing 
of S FR maps, or the adopted conversion factor (Section 
13.21) . Using an ex panded versi o n of t he Monte Carlo mod- 
eling approach of lBlanc et al.l (|2009D , we carry out power 
law fits to our data, avoiding some of the systematic 
biases present in previous work (Section 13. 3|) . Finally, 
we compare our results to a wide collection of literature 
data, demonstrating an emerging consensus with regard 
to the region of parameter space occupied by EsFR-Smoi 
data, if not the interpretation (Section l3.4|) . We conclude 
that in the disks of normal, massive star-forming galax- 
ies, to first order the relationship between Esfr and S„ioi 
can be described by a single depletion time with a factor 
of two scatter. This expands and rein forces the results 
of LOS, BOS, and ISchruba eFall (l20lTI) . 

In Section |4l we show important second-order devia- 
tions from this simple picture, which can easily be missed 
by comparing only Esfr and Smoi- We find systematic 
variations in the molecular gas depletion time, r^"' as 
a function of global galaxy properties (Section 14. 1[) and 
local conditions (Section 14. 2|) . This analysis includes the 
first resolved comparison of r^""' to the local dust-to- 
gas ratio, which is expected to play a key role in setting 
the CO-to-H2 conversion factor, aco and we discuss the 
possibility that aco drives some of the observed r^^p 
variations. We explicitly consider the central regions of 
our targets (Section 14. 4p and show strong evidence for 
lower T^'^, i.e., more efficient star formation, in galaxy 
centers compared to galaxy disks - a phenomenon that 
we discuss in context of recent propo sals for "disk" and 
"starburst" modes of star formation (jPaddi et al.ll201(il : 
IGenzel et al.l [201(11 ). In our best-resolved targets we ex- 
amine how the scatter in t^'^ depends on spatial scale 
(Section 14. 3p . The relationship appears much shallower 
than one would expect for uncorrelated averaging in a 
disk. This suggests either a high degree of large-scale 
synchronization in the star formation process or, more 
likely, widespread systematic, but subtle, variations in 
'^dop "^^6 still-undiagnosed drivers. 

Thus our conclusions may be abstracted to: molecu- 
lar gas and star formation exhibit a first order one-to- 
one scaling but we observe important second order vari- 
ations about this scaling. These include likely conversion 
factor effects, efficient nuclear starbursts, and weak sys- 
tematic variations in t^°' that emerge considering scale- 
dependent scatter or global galaxy properties. The re- 
mainder of this section presents a brief background. Sec- 
tion [2] describes our data and physical parameter estima- 
tion, Sections [3] and |4] motivate these conclusions, and 
Section [5] synthesizes these results and identifies several 
key future directions. 



1.1. Background 

Following iSchmidtl (|1959l . 119631 ). astronomers have 
studied the relationship between gas and the SFR for 
more than 50 years. Most recent work follows Kcnnicut3 
(fT989lll99SbD and compares the surface densities of SFR, 
5^SFR, and neutral (H I-I-H2) gas mass, Egas. Recent work 
focuses heavily on the power law relationship between 
these surface densities (the "Schmidt-Kennicutt law" or 
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the "star formation law"); 



SsFR = A X Eg 



N 



(1) 



An alternative approach tre ats the ratio of gas and 
SFR as the quantity of interest ()Young et al.lll986lll996l 
LOS). This ratio can be phrased as a gas depletion time 
(Ejiioi/Ssfr) or its inverse, a star formation efficiency 
(EgpR/Enioi)- These share the same physical meaning, 
which is the SFR per unit gas. Both convolve a timescale 
with a true efficiency, for example the lifetime of a molec- 
ular cloud with the fraction of gas converted to stars over 
this time. We focus exclusively on molecular gas in this 
paper and phrase this ratio as the molecular gas deple- 
tion time. 



mol 
'dcp 



^SFR 



(2) 



which is the time for star formation to consume the cur- 
rent molecular gas supply. 

The sta t e of t he field is roughly the following. 
lKennicut3 (|199Sb( ) demonstrated a tight, non-linear 
(TV = 1.4 ± 0.15) scaling between galaxy-averaged Esfr 
and Egas spanning from normal disk galaxies to star- 
bursts. Including H I improved the agreement between 
disks and starbursts, but most of the dynamic range and 
the nonlinear slope was driven by the contrast between 
the disks and merger-induced starbursts, especially the 
local ultraluminous infrared galaxy (ULIRG) population. 
The adoption of a single conversion factor for all systems 
also had significant impact; adopting the "ULIRG" con - 
version factor suggested by iDownes fc SolomonI ()1998f ) 
for the lKennicuttr()199Sb() starburst data drives the im- 
plied slope to TV 1.7. Subsequent studies resolved 
galaxy disks — often as radial profiles — and usu- 
ally revealed distinct relationships between EgpR, Ehi, 
and Emoi with shallower indices N f or H2 than H I 
JWong & Blitd 120021: iHever et al.|[200l iKennicutt eFall 
[20071: iSchruba et aLlboill BOS. L08V This suggests that 
the immediate link is between SFR and II2. A more 
aggressive conclusion, motivated by the steep, relatively 
weak relation between Esfr and Ehi is that star for- 
mation in galaxy disks may be broken into two parts: 
(1) the formation of stars in mole cular clouds and (2) 
the balance between H2 and H I (jWong fc Blitd l200a 
iBhtz fc Rosolo wskv 2006. LOS, BOS). 

The fraction of dens e molecular gas also a ppears to 
be a key parameter. iGao fc SolomonI (|2004J ) found a 
roughly fixed ratio of SFR to HCN emission, a dense 
gas tracer, extending from spiral galaxies to starbursts. 
Over the same range, the ratio of SFR to total H2, 
traced by CO emission, varies significantly (though the 
two relate roughly linearly in their normal galaxy sam- 
ple). Galactic studies also highlight the impa.ct of 
density on the SFR on cloud scales (iWu et al.l 120051 : 
IHeiderman etalllMoHLada et al.llMot) . It remains un- 
clear how the dense gas fraction varies inside galaxy 
disks, but me rger-in duced starbursts do show high HCN- 
to-CO ratios (jCao fc Solomon 2004; Garcfa-Buril lo et al.l 

[Ml). 

SFR tr acers and CO can both be observed at high 
redshift. IGenzel et al.l ()2010l ) demonstrated broad con- 
sistency between local H2-SFR relations and those at 
z ~ 1 — 3. One key difference at high-z is the exis- 



tence of disk galaxies with t^I'^ only slightly lower than 
those in local disk galaxies but H2 surface density, Emoi, 
as high as that found in starbursts in the local universe 
(Daddi et al. 2010; Tacconi ct al. 2010;). Merger-driven 
starbursts with similar Emoi can have much lower t^'^, 
suggesting the relevance of another parameter to set t^°' . 
Density and the dynamica l timescale are both good can- 
didates (|Daddi et al.ll20ToHGmi"el et al.ll2010[) . 

Meanwhile, investigations of the Milky Way and the 
nearest galaxies have attempted to connect observed 
scaling relations to the properties of individual star- 
forming regions. These are able to recover the galaxy- 
scale relations at large scale but find enormous scatter in 
the ratios of SFR tracers to molecular ga s on small scales 
(ISchruba et al.l 120101: iChen et al.l 120101: lOnodera etahl 
12010( 1. Detailed studies of Milky Way and LMC clouds 
suggest the time-evolution of individual s tar-forming re- 
gions as a likely sour ce of this scatter fM urravl 120101 : 
iKawamura et al.l |2009() . with th e volume density of in- 
dividual clouds a key parameter (jHeiderman et al.ll2010l : 
iLada et al.llMot ). 

For additional background we refe r the reader to the 
recent review by IKennicutt fc Evan^ ()2012l ) . 

2. DATA 

2.1. Data Sets 

We use HERACLES C0(2-l) maps to infer the distri- 
bution of H2 and GALEX far- ultraviolet (FUV), Spitzer 
infrared (IR) , and literature Ha data to trace recent star 
formation. We supplement these with H I data used to 
mask the CO, derive kinematics, and measure the dust- 
to-gas ratio and with near-IR data used to estimate the 
stellar surface density, E*. 

HERACLES CO: The HERA CO Line Extragalactic 
Survey (HE RACLES) used the H eterodyne Receiver Ar- 
ray (HERA. [Schuster et al.ll2004) on the IRAM 30m tele- 
scope to map CO (2-1) emission from 48 nearby galax- 
ies, of which we use 30 in this paper (see Section [273)) . 
HERACLES combines an IRAM Large Program and sev- 
eral si ngle-semeste r proje cts that spanned from 2007 to 
201 0. iLerov et al] (I2009D presented the first maps (see 
also lSchuster et al]|2007D . The additional data here were 
observed and reduced in a similar manner. The largest 
change is a revised estimate of the main beam efficiency, 
lowering observed intensities by « 10%. This propagates 
to a revised CO (2-l)/(l-0) line ratio estimate, so our es- 
timates of E,noi are largely unaffected compared to BOS 
and LOS. The HERACLES cubes cover out to radii of 
r25 with angular resolution 13" and typical \a sensitivity 
20 mK per 5 km s~^ channel. 

We integrate each cube along the velocity axis to pro- 
duce maps of the integrated intensity along each line of 
sight. We wish to avoid including unnecessary noise in 
this integral and so restrict the velocity range over which 
we integrate to be as small as possible while still con- 
taining the CO line, i.e., we "mask" the cubes. To be 
included in the mask a pixel must meet one of two con- 
ditions: 1) lie within ±25 km s ~^ of the local mea n H I 
velocity (derived from THINGS. [Walter et al.l[200l sup- 
plemented by new and archival H I) or 2) lie in part of 
the spectrum near either two consecutive channels with 
SNR above 4 or three consecutive channels with SNR 
above 3. Condition (2) corresponds to traditional ra- 
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TABLE 1 

Sample 



Galaxy 


D 


res. 


i 


PA 


^25 


'r25 


Multi- 




[Mpc] 


[kpc] 


[°] 


[°] 


['] 


[kpc] 


scale 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 0337 


19.3^ 


1.24"^ 


51 


90 


1.5 


10.6 




NGC 0628 


7.2K 


0.46 


7 


20 


4.9 


10.4 


/ 


NGC 0925 


9.1K 


0.59 


66 


287 


5.4 


14.3 




NGC 2403 


3.2W 


0.21 


63 


124 


7.9 


7.4 


/ 


NGC 2841 


14.1K 


0.91 


74 


153 


3.5 


14.2 




NGC 2903 


8.9W 


0.57 


65 


204 


5.9 


15.2 


/ 


NGC 2976 


3.6^ 


0.23 


65 


335 


3.6 


3.8 




NGC 3049 


19. 2^ 


1.24=* 


58 


28 


1.0 


2.7 




NGC 3184 


11.8K 


0.76 


16 


179 


3.7 






NGC 3198 


I4.1K 


0.91 


72 


215 


3.2 


13.0 




NGC 3351 


9.3K 


0.60 


41 


192 


3.6 


10.6 




NGC 3521 


11. 2K 


0.72 


73 


340 


4.2 


12.9 




NGC 3627 


9.4K 


0.61 


62 


173 


5.1 


13.8 


/ 


NGC 3938 


17.9^ 


1.15'' 


14 


15 


1.8 


6.3 




NGC 4214 


2.9W 


0.19 


44 


65 


3.4 


2.9 




NGC 4254 


I4.4K 


0.93 


32 


55 


2.5 


14.6 




NGC 4321 


14.3^ 


0.92 


30 


153 


3.0 


12.5 




NGC 4536 


14.5^ 


0.94 


59 


299 


3.5 


14.9 




NGC 4559 


7.0K 


0.45 


65 


328 


5.2 


10.7 




NGC 4569 


9.86^ 


0.64 


66 


23 


4.6 


26.5 




NGC 4579 


16. 4K 


1.06'' 


39 


100 


2.5 


15.0 




NGC 4625 


9.3K 


0.60 


47 


330 


0.7 


1.9 




NGC 4725 


II.9K 


0.77 


54 


36 


4.9 


13.2 


/ 


NGC 4736 


4.7K 


0.30 


41 


296 


3.9 


5.3 


/ 


NGC 5055 


7.9K 


0.51 


59 


102 


5.9 


17.3 




NGC 5194 


7.9W 


0.52 


20 


172 


3.9 


9.0 


/ 


NGC 5457 


6.7K 


0.43 


18 


39 


12.0 


25.8 


/ 


NGC 5713 


21. 4K 


1.38=' 


48 


11 


1.2 


9.5 




NGC 6946 


6.8K 


0.44 


33 


243 


5.7 


9.8 


/ 


NGC 7331 


14.5^ 


0.94 


76 


168 


4.6 


19.5 





Note. — Sample used in this paper. Columns give (1) galaxy 
name; (2) adopted distance in Mpc; (3) FWHM spatial resolu- 
tion of HERACLES data at that distance, in kiloparsecs; (4) 
adopted inclination and (5) position angle in degrees; adopted 
radius of the the B-band 25"' magnitude isophote, used to nor- 
malize the radius in (6) arcminutes and (7) kiloparsecs. Most 
analysis in this paper considers data inside 0.75 r25. Column 
(8) indicates if the galaxy is cl ose and large enough for the mul- 
tiscale analysis in Section 14.31 

" Too distant to convolve to 1 kpc resolution. Included in anal- 
ysis at native resolution. 

K,W Distanc e adop ted from K: IKennicutt et aO l|20111 1 or W: 
IWalter et all ([2001). 

dio masking fe.g- IHelfer eraI|[2003HWalter et al.|[2008l ). 

Condition (1) is less conventional, but important to our 
analysis. Integrating over the H I line, which is detected 
throughout our targets, guarantees that we have an inte- 
grated intensity measurement along each line of sight, 
even lines of sight that lack bright CO emission (see 
iSchruba et al.l 1201 IL for detailed discussion of this ap- 
proach). This avoids a traditional weakness of masking, 
that nondetections are difficult to deal with quantita- 
tively. We calculate maps of the statistical uncertainty 
in the integrated CO intensity from the combination of 
the mask and estimates of the noise derived from signal- 
free regions. The result is an integrated intensity and 
associated uncertainty for each line of sight in the HER- 
ACLES mask. 

SINGS and LVL IR: We use maps of IR emission from 
3.6-160/im from the Spitzer Infrared N earby Galaxies 
burvey (SINGS, IKennicutt et al.l l2003aD and the Local 



Volume Legacy survey (LVL, iDale et al.ll2009D. W e de- 
scribe the processing of these maps in lLerov etalllfMl . 
hereafter L12). 

SINGS, LVL, and Literature Ha: Both SINGS and 
LVL published continuum-subtracted Ha images for 
most of our sample. We supplement these with literature 
maps, particularly from the GoldMine and Palomar-Las 
Campanas surveys. L12 describe our approach to these 
maps (masking, N II correction, flux scaling, background 
subtraction) and list the source of the Ha data for each 
galaxy. 

GALEX UV: For 24 galaxies, we use NUV and 
FUV maps from the Nearby Galaxy Survey (NGS, 
IGil de Paz erall[200l . For one galaxy, we use a map 
from the Medium Imaging Survey (MIS) and we take 
maps for five targets from the All-sky Imaging Survey 
(AIS). L12 describe our processing. 

BIMA+12-m and NRG 45-m GO (1-0) Maps: A sub- 
set of o ur targets ha.ve als o been observed by the BIMA 
SONG (IHelfer et al.ll2003D or the Nobevam a CO Atlas of 
Nearby Spiral Galaxies (|Kuno et al.ll2007 [). Where these 
data are available, we apply our HERACLES masks to 
these maps and measure CO (1-0) intensity. Wc only use 
BIMA SONG maps that include short-spacing data from 
the Kitt Peak 12-m. 

THINGS and Supplemental Hi.- We assemble H I maps 
for all targets, which we use to mask the CO, estimate the 
dust-to-gas ratio, explore 24/im cirrus corrections, and 
derive ap proximate rot ation curves. These come from 
THINGS ([Walter et aII[200Sl and a collection of new and 
archival VLA data (programs AL731 and AL735). These 
supplemental H I are C-l-D configuration maps with res- 
olutions 13"-25". We reduced and imaged these in a 
standard way using the CASA package. 

2.2. Physical Parameter Estimates 

Following standard practice in this field, we estimate 
physical parameters from observables. Despite the in- 
trinsic uncertainty involved in this process, these esti- 
mates play a fundamental role in enhancing our under- 
standing of the physics of galaxy and star formation, 
as demonstrated from the earliest works in this subfield 
lYoung & Knczck 1989; Kennicutt 1989). We adopt an 
approach largely oriented to physical quantities, but dis- 
cuss the impact of our assumptions throughout. 

CO Intensity to H2: We convert CO (2-1) intensity to 
H2 mass via 

S]^oi[Mqpc-2]^6.3 (^-^) (^^) /CO [Kkms-i], 

(3) 

where R21 is the CO(2-l)-to-CO(l-0) line ratio and 
aco is the CO(l-0)-to-H2 conversion factor. By de- 
fault, we adopt a Galactic conversion factor, — 

4.35 Mq pc~2 (Kkms"^)-i equivalent to X co = 

2 X 10^°cm-^ (K kms-^)-^ (jStrong fc Mattoxl Il996l : 
IDame et all 120011) and a line ratio of R21 = 0.7. This 
line ratio is slightl y lower than the R21 =0.8 derived by 
iLerov etldl ()2009[ ). reflecting the revised efficiency used 
in the reduction. The appendix motivates this value us- 
ing integrated flux ratios and follow-up spectroscopy of 
HERACLES targets. Equation [3] and aU "Emoi" in this 
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paper include a factor of 1.36 to account for helium. Be- 
cause we consider only molecular gas, any results that 
we derive using a fixed aco can be straightforwardly re- 
stated in terms of CO intensity. 

We adopt this "Galactic" to facilitate clean com- 
parison to previ ous work, but improved estimates exist 
for HERACLES. ISandstrom eTalT (|20ll) solved directly 
for the CO-to-H2 conversion factor across the HERA- 
CLES sample using dust as an independent tracer of 
the gas mass. They find a somewhat lower average 
•^co" * We quote this as a CO (1-0) conversion 

factor, though ISandstrom et al.l ()2012l ) directly solve for 
the CO (2-1) conversion factor. They find a C0(2-l) 
conversion factor of a^o^ w 4.4 Mq pc-2 (Kkms"^)^i, 
compared to our "Galactic" C0(2-l) conversion factor 

(fro = 6.3 Mq pc-2 (Kkms"V^- ISandstrom et al.l 
(|2012l) find « 0.4 dex point-to-point scatter, of which 
« 0.3 may be intrinsic with the remainder solution uncer- 
tainties. Because [Sandstrom et al. (2012) solve directly 
for a CO (2-1) conversion factor using the same HERA- 
CLES data employed in this paper, these values should 
be borne in mind when reading our results. Our results 
remain pinned to a Galactic CO(l-O) conversion factor of 
a\r^ = 4.35 M0 pc-2 (Kkms"^)-i that may be w 30% 
too high, on average. As a result, a systematic bias of 
w 30% in Smoi appears plausible with factor of two vari- 
ations in the conversion factor point-to-point. 

In addition to a fixed conversion aco j we consider the 
effects of variations in aco due to decreased dust shield- 
ing at low metallicity and variations in the linewidth, 
optical depth, and temperature of CO in galaxy cen- 
ters. Our "va ri able" aco builds on the work of 
ISandstrom et al.l (|2012t ). who compare HI, C0(2-l) 
and E n„st in 22 HERACLES galaxies and iWolfire et all 
([2010') , who consider the effects of dust shielding on the 
"CO-dark" layer of molecular clouds, where most H is 
H2. The aco prescription combines three terms 

aco = "CO CCO-dark (£'/G) Ccontcr (''gal) ■ (4) 

Here a^o — 6-3 M© pc^^ (Kkms^^)^^ is our fiducial 
CO (2-1) conversion factor in the disk of a galaxy at 
solar metallicity (Equation [3]). The term cco-dark rep- 
resents a correction to the H2 mass to refiect the H2 
in a CO-dark layer not directly tr aced by CO emission . 
We calculate this factor following IWolfire et al.l ()2010( ). 
assuming that all GMCs share a fixed surface density, 
(Egmc) and adopting a linear scalin g between the dust- 
to-ga s ratio and metallicity (see also lGlover fc Mac Lowl 
|2011|) . In this case 

cco-dark {D/G') ^ 0.65 exp ■ (5) 

Here D/G' is normalized to our adopted "Galac- 
tic" value of 0.01, with the normalization constructed 
to yield aco = Q^co for D/G" = 1. Sioo = 
(Sgmc) /lOO M0 pc~^. The appendix presents this cal- 
culation in detail. 

We consider two cases: (Sgmc) — 100 M© pc^^ 
("E = 100") and (Sgmc) = 50 M© pc'^ ("S = 50"). 
"E = 100" refiects a typical surface density that is of- 
ten assumed and observed for extragalactic GMCs (e.g.. 



iBolatto et all 120081: iNaravanan et all I2012D . O ver the 

range of D/G' that we consider, this prescription rea- 
sonably resembles the shallow power law dependences 
of aco on metallicity cal culated from siinulatio ns by 
iFeldmann et al.l (|2012l i and lNaravanan et al.l (|2012| ) , who 
both suggest aco ~ "(Sgmc) = 50 M© pc'^" 

yields a steeper dependence of aco on metallicity over 
our range of interest but we will see in Sections 14.1! and 
14.2! that it offers a simple way to account for most of the 
dependence of on dust-to-gas ratios in our observa- 
tions. Low surface density GMCs or significant contribu- 
tion of "translucent" [Ay ^ 1-2 mag) gas to the overall 
CO emission ar e supported by observations of the Milky 
Wav lHever eFa l. (2009) and Liszt et al. (2010), L MC ob - 
servations by " Hughes et all ([20101); iWong et all (|201lD . 
and M31 (Schruba et al. in prep.) but may not be ap- 
propriate for more actively star-forming systems (Hughes 
et al., submitted). We return to this issue in Section [H 

We calculate aco{D/G') using D/G' derived for fixed 
aco- Given observations of Ehi, Ico, Sdust, and a pre- 
scription for acoiD/G'), one can simultaneously solve 
for aco and D/G'. The solution is often multivalued 
and unstable, though not intractable. However, after ex- 
perimentation and compa rison with the se lf -consi stent 
Herschel-hased results of ISandstrom et al.l (|2012D . we 
found that the process does not clearly improve our esti- 
mates. In the interests of clarity and simplicity, we work 
with only D/G calculated using fixed aco throughout 
the paper. This simplification biases our aco estimate 
high by « 8% ("E = 100") and « 15% ("E = 50"). 
With improved Edust estimates, we expect that the self- 
consistent treatment will become necessary. 

The third term, Cccntor, accounts f or depressed val- 
ues of aco in the centers of galaxies. ISandstrom et al.l 
(|2012f l find such depressions in the centers of many 
systems (see also llsraeli l2009al lbl). These likely re- 
flect the same line-broadening and temperature effects 
that dri ve the commonly in voked "ULIRG conversion 
factor" ([D ownes fc Solomon 1998), though the depres- 
sion observed by ISandstrom et al. (2012) have lower 
magnitude than th e factor of 5 depression fo u nd by 
iDownes fc Solomon! (|1998l ). ISandstrom et al.l ()2012[ 1 
could not identify a unique observational driver for 
these depressions, though they correlate well with stel- 
lar surface density. Instead, they appear to be present 
with varying magnitudes in the centers of most systems 
with bright central CO emission. Following their rec- 
ommendation, we apply this correction where Tgai < 
0.1r25 in systems that have such central CO concentra- 
tions. Whenever availab le, we adopt Cccntor directly from 
ISandstrom et al.l ([2012), taking the factor by which the 
central aco falls below the mean for the disk of that 
galaxy. For systems with central CO c o ncent rations but 
not in the sample of ISandstrom et aTl (|2012l ). we apply 
a factor of two depression, again following their recom- 
mendations. The appendix presents additional details. 

SFR from Ha, UV, and 24^^ni Emission: L12 com- 
bined UV, Ha, and IR emission to estimate the recent 
star formation rate surface density, Esfr, at 1 kpc res- 
olution (the limiting common physical resolution of the 
HERACLES survey) for our sample. We adopt their es- 
timates and refer the reader to that work for detailed 
discussion. Briefly, our baseline estimate of Esfr com- 
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bines Ha and infrared emission at via 



on comparison to iZibetti et al.l ()2009D , we use 



SsFR [Mq yr-1 kpc-2] =634 /Ha [erg s sr^^] + (6) 

0.00 3 25 /24 [MJy si~\}) 

where /hq and /24/im refer to the line-integrated Ha in- 
tensity and intensity at 24/im. 

The Ha emission captures direct emission from H II 
regions powered by massive young stars while the 24/im 
emission accounts for recent star formation obscured by 
dust. Before estimating Ssfr, we correct our 24/im maps 
for the effects of heating of dust by a weak, pervasive 
radiation field (i.e., a "cirrus") with magnitude derived 
from modeling the infrared spectral energy distribution. 
The cirrus removed corresponds to the expected emis- 
sion from the local dust mass illuminated by a quiescent 
radiation field, typically ~ 0.6 times the Solar neigh- 
borhood interstellar radiation field (see L12 for details). 
We derive the appropriate weighting for the combination 
of Ha and 24/im emission based on comparing our pro- 
cessed Ha and 24/xm maps to literature estimates of Ha 
extincti on. The resulting line ar combination resembles 
that of iKennicutt et all (|2007D but places slightly more 
weight on the 24/zm term. For comparison, we also esti- 
mate EsFR from combining FUV and 24/xm emission and 
taking Ha alone while assuming a typical 1 magnitude 
of extinction. 

L12 estimate a substantial uncertainty in the absolute 
calibration of "hybrid" UV-flR or Ha-I-IR tracers, with 
magnitude » 50%. In addition to this overall uncertainty 
in the calibration, they derive a point-to-point uncer- 
tainty in SsFR of ~ 0.15 dex based on intercomparison 
of different estimates. 

Dust Properties: In order to measure dust properties, 
we convolve the Spitzer 24, 70, and 160/im data and the 
CO and H I maps to the resolution of the Spitzer 160/xm 
data. At this resolution, we build ra dial profiles of each 
band and then fit the dust models of lDraine fc Lil (|2007[ ) 
to these profiles. These fits, presented in L12, provide 
us with radial estimates of the dust-to-gas ratio, D/G, 
and are used to help account for "cirrus" contamination 
when estimating Ssfr- Note that the ~ 40" resolution 
of the 160/im data used to measure these dust proper- 
ties is significantly coarser than the 1 kpc resolution used 
for the rest of our data. Where possible, we have com- 
pared our Spitzer-hased dust masses to masses estimated 
using the improved SEP coverage offered by Herschel 
(;e.g.. lAniano et al.|[20T2l ): above Edust « 0.05 Mq pc'^ 
the median offset between the Herschel and Spitzer based 
dust masses is only « 10%; however the dust masses de- 
rived for individual rings using only Spitzer do scatter 
by ~ 0.3 dex (a factor of two) compared to Herschel- 
based dust masses and show weak systematic trends with 
the sense that Spitzer underestimates the mass of cooler 
(~ 15 K) dust in the ou tskirts of galaxies (both consis- 
tent with the analysis of lDraine et al.|[2007f) . We expect 
that once Herschel images become available, they will 
significantly improve the accuracy of dust-based portion 
of this analysis. 

Stellar Mass: To estimate the stell ar mass for whole 
galax ies, we draw 3.6/im fluxes from iDale et all (120071 
120091) ■ convert to a luminosity using our adopted dis- 
tance, and apply a fixed 3.6/xm mass-to- light ratio. Based 



[Mo pc-2] = 200 /3.6 [MJy sr'^] 



(8) 



which is ~ 30% lower than LOS. This value is uncertain 
by - 50%. 

We estimate the stellar surface density, E*, for each 
kpc-sized element from the contaminant-corrected 3.6/xm 
maps of M eidt et"al] ()2012D . Starting from a reprocessing 
of the SINGS data (as part of the S4G survev lSheth et al.! 
I2010f) . they used independent component analysis to re- 
move contamination by young stars and hot dust from 
the overall maps. These contaminants make a minor 
contribution to the overall 3.6/^m fiux but may be im- 
portant locally. We convert the contaminant-corrected 
3.6/im maps to E* estimates using Equati on [5] 

Rotation Velocities: Following LOS and IBoissier et al.l 
(2003) we work with a simple two-parameter fit to the 
rotation curve of each galaxy 



Vrot (?'gal) = Vfla,t 



1 — exp 



flat 



(9) 



with Wfiat and Zfiat free paramete rs. We deri ve these from 
fits to the rotation curves of de Blok et all ^008) wher- 
ever they are available. Where these are not available, 
we carry out our own tilted ring fits to the combined H I 
and CO first moment maps. We use these fits to calcu- 
late the orbital time Torb = 27rrgai/i'rot for each line of 
sight. 

2.3. Sample and Galaxy Properties 

We present measurements for galaxies meeting the fol- 
lowing criteria: 1) a HERACLES CO map containing a 
clear CO detection (S/N> 5 over a significant area and 
multiple channels), 2) Spitzer data at 24/im, and 3) in- 
clination < 75°. The first condition excludes low mass 
galaxies without CO detections (these are discussed in 
iSchruba et al.|[2012[ ). The second removes a few targets 
with saturated or incomplete Spitzer coverage. The third 
removes a handful of edge-on galaxies. We are left with 
the 30 disk galaxies listed in Table [TJ 

For each target, Table [T] gives the distance, phys- 
ical resolution of the HERACLES maps at that dis- 
tance, inclination, position angle, and optical radius. 
The Table notes note the subset of galaxies that that 
are close and large enough for us to carry out the 
multi-resolution analysis in Se c tion |4.3I We adopt 
distances from IKennicutt et al.l ()201ll ) where possible 
and from iWalter et al.l (120081) elsew here. We take ori- 
entations from 'Waher et al.' (200S') and from LEDA 
([Prugnici & Hcraudcau 199S) and NED elsewhere. 

Table [2] reports integrated and disk-average properties 
for our sample. We report our integrated stellar mass es- 
timate, galaxy morphology, metallicity and dust-to-gas 
ratio at « 0.4 r25, average gas mass and star forma- 
tion rate surface density inside 0.75 r25, our parameter- 
ized rotation curve fit, an d the orbital time at 0.4 r25. 
We take metallicities from IMoustakas et al.| (j2010D , av- 
eraging their PT05 and KK04 strong-line calibrations. 
They argue that these two calibrations bracket the true 
metallicity and that t he relative ordering of metallici- 
ties is robust (see also lKewlev fc Ellisonll2008l) . but the 
uncertainty in the absolute value is considerable. For 
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cases where iMoustakas et al.l (|201ClD do not present a 
me tallicity, we draw o ne fr om the recent comp ilations 
bv lMarble et~all (|2010l ) and lCalzetti"eFaII (|20lol ). 

2.4. Methodology 

We sample our targets at 1 kpc resolution. This is fine 
enough to isolate many key physical conditions in the in- 
terstellar medium (ISM): metallicity, coarse kinematics, 
gas and stellar surface density. At the same time, we 
expect to av erage several star for ming regions in each el- 
ement (e.g.. lSchruba et al.|[20Tot) . with M^oi > lO'^ Mq 
and M* > 10^ Mq formed over the last ~ 5 Myr in each 
element. This minimizes concerns about evolution of in- 
dividual regions, sampling the IMF, and drift of stars or 
leakage of ionizing photons from their parent region. 

We convolve each map to have a symmetric gaussian 
beam with FWHM 1 kpc. For the Spitzer 24/im maps 
we first convert from the MIPS PSF to a 13" gaussian 
beam using a kernel kindly provided by K. Gordon, then 
we convolve to 1 kpc. This exercise effectively places 
our targets at a common distance but does not account 
for foreshortening along the minor axis. Five galaxies 
are too distant to convolve to 1 kpc. We mark these in 
Table [T] and include them in our analysis at their native 
resolution. 

We sample each map to generate a set of intensity mea- 
surements. The sampling points are distributed on a 
hexagonal grid with points spaced by 0.5 kpc, one half- 
resolution element. At each sampling point we measure 
C0(2-l) intensity. Hi intensity, a suite of star formation 
rate tracers (described in L12), dust properties, and S,. 
We use these to estimate physical conditions as described 
above and in L12, taking into account the inclination of 
the galaxy. 

We also identify a sample of galaxies to study the ef- 
fects of physical resolution. Nine galaxies, marked in 
Table [U have both the proximity and extent to allow us 
to test the effect of physical resolution on our results. 
We convolve these to a succession of physical resolutions 
from 0.6 to 2.4 kpc for further analysis (Section 14. 3p . 

We treat regions with Esfr < 10"'^ Mq yr^^ kpc^^ or 
Ico < 2.5 X (Tco as upper limits and consider only points 
with Tgai < 0.75 r25 — t he HERACLES maps contain sig- 
nal outside this radius (iSchruba et al.ll2Qlll ) but mostly 
not significant emission over individual lines-of-sight. In 
total we have ^ 14, 500 lines of sight with at least one 
significant measurement, of which 1, 900 have CO upper 
limits and 1,650 have SF upper limits. Points for which 
neither measurement is significant are not considered in 
the analysis. Nyquist sampling the maps in a hexago- 
nal pattern leads to oversampling by a factor of ~ 5, 
so that this corresponds to > 2, 000 independent mea- 
surements. The maximum (2.5a) upper limit on Smoi is 
w 6 M0 pc~^, the median upper limit is w 2.6 M0 pc~^. 

2.5. Literature Data 

We compare our results to recent measurements of SFR 
and molecular gas. These employ a variety of sampling 
schemes and SFR tracers. We adjust each to match our 
adopted CO-to-H2 conversion factor and IMF. Contrast- 
ing our approach with these data illuminates the impact 
of methodology and allows us to explore whether diverse 
observations yield consistent results under matched as- 



sumptions. 

iKennicuttI (|1998bD presented disk-averaged measure- 
ments for 57 normal spiral galaxies and 15 starburst 
galaxies. He used literature CO with a fixed aco to 
estimate Smoi- To estimate Ssfr, he used Ha in disk 
g alaxies and IR emissi on in starbursts. 

iCalzetti et al.l ()2010[ ) estimated disk-averaged Ssfr for 
a large set of nearby galaxies. We cross-inde x these 
with integra t ed C O flux es from lYoung et al.l ([T995), 
IHelfer et all ()2003[ ). and ILerov et al.l (|2009D to derive 
Einoi assuming that CO emission covers the same area 
as Ha. From the combination of these data we have 
disk- average Ssfr a nd Sm oi estimates f or 41 galaxies. 

Sai ntonge et al.l (|2012D . following iSaintonge et al.l 
(|2011i) . present the COLDGASS survey, which obtained 
integrated molecular gas mass and SFRs for 366 galax- 
ies with > 10^° M0, 215 with secure CO detections. 
This large survey represents the best sample of integrated 
galaxy measurements to date. To convert to surface den- 
sities, we take the area o f the star-forming disk in these 
galaxies to be 0.75 r25. ISaintonge et al.l (|2012D derive 
their SFRs from SED modeling that yields results close to 
what one would obtain converting the UV-I-IR luminos- 
ity directly to a SFR. This yields higher SFRs than our 
approach for matched measurements. Comparing galax- 
ies with matched stellar mass or molecular gas content, 
we find the offset to be « 0.19 dex, a factor of w 1.55. 
This agrees well with what one would expect account- 
ing for our subtraction of an IR cirrus with magnitude 
« 1.2 and our 24/im coefficient, which is « 1.2 lower than 
what one would adopt to match a bolometric TIR SFR 
i ndicator (see L1 2 for calculations and discussion). 

ILerov et"aLl (I2005D combin ed new data with 
measure ments bv | You ng "eFaT (T995^, 'Elfhag^eialJ 
(il99d). iTav l or et al.l (11998). Bokcr et al. (2003), and 
IMurgia et al.l (12002) to compare Ssfr and Emoi for 
individual ~ 30-50" pointings in a wide sample of 
nearby galaxies. They estimat e Ssfr from the 20cm 
radio continuum (jCondonl 119921 ). These low- resolution 
pointings typically cover several kpc, a larger area than 
our resolution elements but less than an average over a 
whole galaxy disk. 

Wong & Blitz (20021) iSchuster et al.l (|2007[) . and 
ICrosthwaite fc Turneii ()2007D presented radial pro- 
files of Smni and S rfr. for several nearby galaxies. 
IWong fc Blitd (j2002 ^ target ed 7 nearb y spirals, using 
Ha to calculate Esfr- S chuster et al.l (|2b07.1 targeted 
M51 and derived Esfr from 20-cm ra dio continuum 
to est imat e Esfr - We only pr esent the iWong fc Blitd 
(|2002D and lSchuster et al.l (|2007D profiles down to Emoi ~ 
5 Mq pc~^, below which we con sider them some- 
what unreliable. ICro sthwa ite fc Turne r (2007) targeted 
N GC 6946 and use d I R emi ssion to estimate Esfr- 

iKennicutt et all |2007D . iBlanc e t al.l (I2009I) . 
I Rahman et al.l pOllf T and [Rahman et al.l ()2012f ) 
targeted small re gions, similar to BOS a nd the work 
presented here. IK ennicutt et al.l (|2007f ) focused on 
luminous regions in M51, mainly in the spirals arms. 
They infer Esfr from a c o mbina tion of Ha and 24 /xm 
emission. iRahman et al.l (|2011[ ) explored a range of 
methodologies. We focus on their most robust measure- 
ments, drawn from bright regions in NGC 4254 with 
Esfr, from a combination of NUV and 24 /im emission. 
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TABLE 2 

Sample Properties 



Galaxy 


log(M.) 


Morphology 


z 


D/G 


(SHI+H2) 


(SsFR> 


■"flat 


'flat 


(forb) 




loffin \'Mr^^ 


T-TvDe 


fl2+lofffO/Hll 






[10-3 M0yr-i 


[km s^^] 


fkDcl 


flO* vrl 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


JNCjO Vooi 


9.9 


6.7 


8.51 


0.004 


21 


14 


130 


2.9 


2.4 




10.0 


5.2 


8.69 


0.012 


15 


4.0 


200 


0.8 


1.3 




9.7 


7.0 


8.52 


0.004 


7.5 


1.3 


140 


6.9 


4.4 


■\T o /I r\ Q 


9.6 


6.0 


8.57 


0.009 


10 


3.3 


120 


0.95 


1.6 




10.7 


3.0 


8.88 


0.037 


4.6 


1.4 


310 


2.3 


1.2 


■\Ti^/~i onr\o 


10.4 


4.0 


8.90'^ 


0.012 


12 


5.7 


210 


2.4 


2.0 


InCjC zy/o 


9.0 


5.2 


8.67 


0.008 


7.6 


4.4 


88 


1.1 


1.4 


■\T/^/-i on /in 
In <jU oU4y 


9.5 


2.5 


8.82 


0.005 


8.0 


10 


180 


3.0 


1.5 




10.2 


6.0 


8.83 


0.018 


14 


2.8 


200 


2.5 


1.8 


■\T/^ /-I o T no 


10.0 


5.2 


8.62 


0.012 


8.4 


2.3 


150 


3.0 


2.6 


■\T o o c T 


10.1 


3.1 


8.90 


0.018 


8.5 


5.2 


200 


1.1 


1.2 




10.7 


4.0 


8.70 


0.012 


22 


7.8 


229 


1.5 


1.5 




10.5 


3.1 


8.67 


0.016 


13 


7.7 


190 


1.1 


1.8 


■\Ti^/~i onoo 


10.3 


5.1 


8.71^^ 


0.018 


22 


7.9 


140 


0.73 


1.6 


INijL/ 4zi4 


8.7 


9.8 


8.36^ 


U.OU08 


9.2 


8.4 


350 


11 


2.0 


4254 


10.5 


5.2 


8.79 


0.01 


47 


18 


170 


1.4 


1.6 


INCjO 4ci21 


10.6 


4.1 


8.84 


0.012 


30 


9.0 


229 


1.8 


1.4 


NijL> 4ooo 


10.2 


4.3 


8.61 


0.005 


13 


6.8 


180 


0.7 


2.0 


InCtU 4ooy 


9.5 


6.0 


8.55 


0.005 


10 


1.8 


100 


2.1 


2.9 


NGC 4569 


10.2 


2.4 


8.88^^ 


0.017 


8.5 


1.9 


220 


3.2 


1.8 


NGC 4579 


10.7 


2.8 


8.93^= 


0.021 


13 


3.8 


270 


1.7 


1.2 


NGC 4625 


8.9 


8.8 


8.70 


0.011 


8.5 


6.6 


27 


0.53 


2.3 


NGC 4725 


10.5 


2.2 


8.73 


0.03 


4.5 


0.75 


220 


1.1 


1.9 


NGC 4736 


10.2 


2.4 


8.66 


0.008 


17 


10 


170 


0.25 


0.77 


NGC 5055 


10.5 


4.0 


8.77 


0.02 


18 


4.1 


200 


0.71 


1.7 


NGC 5194 


10.5 


4.0 


8.87 


0.02 


53 


20 


210 


0.58 


1.0 


NGC 5457 


10.4 


6.0 


8.46^= 


0.013 


10 


2.4 


210 


1.2 


2.7 


NGC 5713 


10.3 


4.0 


8.64 


0.006 


54 


37 








NGC 6946 


10.5 


5.9 


8.73 


0.007 


37 


21 


190 


1.2 


1.5 


NGC 7331 


10.8 


3.9 


8.68 


0.01 


16 


4.4 


260 


1.9 


1.9 



Note. — Properties of samp le galaxies. Columns give (1) galaxy name; (2) integrated stellar mass o f whole galaxies based o n 
3.6^m flux of I Dale et al.l 12007. 2009); (3) morphology; (4) "characteristic" metallicity at 0.4 r25 from [Moustakas et al.l | |2010I '|. 
averaging their PT05 and KK04 calibrations; (5) dust-to-gas ratio at 0.4 r25 based on our modeling of Spitzer data; (6) average 
H I-(-H2 surface density inside 0.75 r2B', (7) average star formation rate surface density inside 0.75 r25; parameters for simple 
rotation curve fit, (8) iifjat ^^nd (9) (flat; and (10) orbital time at 0.4 r25 based on the rotation curve. 

" MetalUcity from compilation of lCalzetti et al.l | |20TO|'| and lMarble et al.l i20Wl ) or lKennicutt et al.l 12003 (NGC 5457). 



I Rahman et alJ (|2012 fl extended this work to consider the 
Ml set of CARMA STING galaxies, using only the 24^m 
emission wi t h a nonlinear calibration to infer Ssfr- 
iBlanc et all ()2009D studied the central 4.1 x 4.1 kpc^ of 
M51, deriving Esfr from Ha spectroscopy corrected 
using the Balmer decrement. 



3. S; 



SFR-^MOL 



SCALING RELATIONS: FIRST ORDER 
CONSTANCY OF r^OL 



We estimate Ssfr and Smoi for ^ 14, 500 points in 30 
nearby galaxies. In this section we analyze these data in 
the context of a traditional "star formation law" scaling 
relation ( jjl.ip . We show the data distribution in Ssfr- 
Smoi parameter space (^ 13.11) and examine how this de- 
pends on methodology (j }3^. Using a Mon te Carlo tech- 
nique based on that of Bl anc et al.l (j2009D . we consider 
the best fit power-law to the ensemble data and individ- 
ual galaxies (t j3.3|) . We compare our results to a broad 
sample of literature data (! j3.4p . 

3.1. Combined Measurement 

Figure [T] compares Ssfr, estimated from Ha-|-24/im, 
and Smoi at 1 kpc resolution for our whole sample. Indi- 
vidual kpc resolution lines of sight appear as gray points 
and the red points show the median Esfr and standard 
deviation after binning the data by Smoi- In the top right 



panel and bottom row, blue contours show data density 
adopting different weightings. The top right panel gives 
identical weight to each line of sight, treating each kpc^ 
as equal regardless of location. The bottom left panel 
gives equal weight to each galaxy and so weights mea- 
surements from small galaxies with little area more than 
measurements from large galaxies. The bottom right 
panel treats each radial ring in each galaxy equally, and 
so gives more weight to the central parts of galaxies than 
their outer regions. Dashed lines here and throughout 
this paper indicate fixed r^"' and a horizontal line indi- 
cates the limit of our Esfr measurements. In the top left 
panel, dark points show measurements where one quan- 
tity is an upper limit. Table [3] summarizes key values 
from the plots in this section. 

The top rows of Figure [T] and Table [3] show the good 
correspondence between Esfr and Emoi that we have 
previously found in the HERACLES s ample (B08,L08, 
LSchruba et al. 201ll: IBigiel et al.ll20lil) . Our dynamic 
range at 1 kpc resolution spans from Esfr ~ 10~^ 



to IQ- 



kpc and Emoi from a few to 



100 M© pc~ . Across this range, E,„oi and Esfr cor- 
relate well, exhibiting a Spearman rank correlation coef- 
ficient > 0.7 for most tracers and weightings. This quan- 
tifies the tight, one-to-one relationship visible by eye in 
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Fig. 1. — Star formation rate surface density, Espr, estimated from Ha+24/im emission, as a function of molecular gas surface density, 
Smol, derived from CO(2-l) emission for 30 nearby disk galaxies. The top left panel shows individual points (dark gray points show upper 
limits) with the running median and standard deviation indicated by red points and error bars. The red points with error bars from the first 
panel appear in all four panels to allow easy comparison. Dotted lines indicated fixed H2 depletion times; the number indicates logjg TQ^p 
in yr. The top right panel shows the density of the data in the top left panel. In the bottom panels we vary the weighting used to derive 
data density. The bottom left panel gives equal weight to each galaxy. The bottom right panel gives equal weight to each galaxy and each 
radial bin. 



the top row. 

The median t^I'^ weighting each line of sight equally 
is 2.2 Gyr with a scatter of 0.3 dex, a factor of two. 
The absolute value of the median t^^, i.e., the scale 
of the X- and y-axes in Figure [1] depends on the cali- 
bration of our SFR tracer and CO-to-H2 conversion fac- 
tor. Each remains uncertain at the 30-50% level and we 
suggest that an overall uncertainty of 60% on the ab- 
solute value of T^p' represents a realistic, if somewhat 
conservative, value. Our Emoi and Ssfr estimates can 
be compa red internally with mu ch better accuracy than 
this fL12. lSandstrom et al.|[20T^ . so we suggest that this 
uncertainty be viewed as an overall scaling of our results. 

The bottom row in Figure [T] begins to reveal the devi- 
ations from a simple one-to-one scaling that will be the 
subject of Section m Weighting all galaxies equally (bot- 



tom left panel) reveals a significant population of low 
Smoi, high SsFR, low r^°' data. This drives the median 
depletion time for the sample from w 2.2 Gyr, weighting 
by line-of-sight, to w 1.3 Gyr, weighting by galaxy. In 
Section 14.11 we show that these low apparent t^^"' orig- 

inate from low-mass, low-metallicity systems (see also 
Schruba et aT1l20TlHK7^holz et al.ll201ll : ISchruba et al.l 
2012| ). Because of their small size, these systems do not 
contribute many data compared to large, metal-rich spi- 
rals. Therefore, they only weakly influence the overall 
data distribution seen in the top row. We examine t^"^ 
as a function of host galaxy properties and local condi- 
tions in Sections 14. II and 14. 21 In the appendix we present 
^SFR — Smoi relations for individual galaxies (see also Ta- 
ble [2|) , allowing the reader to see how Figure [1] emerges 
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TABLE 3 

T^"°' AT 1 KPC Resolution 



from the superposition of individual systems (see also 
Section 13. 3p . 

Weighting radial rings equally (bottom right panel) 
highlights these same low r^p-low Emoi galaxies and 
brings out an additional low r^"' population at higher 
Smoi- These point emerge because the radial weighting 
emphasizes points in the central parts of galaxies rela- 
tive to their outskirts. We show in Section that the 
central regions of many of our targets exhibit enhanced 
efficiency compare to their disks. As will small galaxies, 
these central regions contribute only a tiny fraction of 
the area in our survey and thus exert little impact on 
the plots in the top row. 

Figure [1] thus illustrates our main conclusions: a first 
order simple linear correlation between Ssfr and Smoi 
and real second order variations. It also illustrates 
the limitation of considering only SspR-Smoi parameter 
space to elicit these second-order variations. Metallic- 
ity, dust-to-gas ratio, and position in a galaxy all play 
key roles but are not encoded in this plot, leading to 
double-valued Ssfr at fixed Smoi in some regimes. We 
explore these systematic variations in r^"' and motivate 
our explanations throughout the rest of the paper. 

3.2. Relationship for Different SFR and Molecular Gas 

Tracers 

Figure [1] shows our best-estimate Ssfr and Smoi com- 
puted from fixed aco- Many approaches exist to esti- 



mate each quantity (see references in lLerov et al.|[201lL 
L12) and the recent literature includes many claims 
about the effect of physical parameter estimation on the 
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relation between Ssfr, and Ej, 



In this section, we 



Note. — Median molecular gas depletion time, scatter, and corre- 
lation between EsFR. a-nd Smo in our sample. Line-of-sight averages 
treat each kpc-resolution line of sight as equal. Galaxy averages re- 
fer to T^°' = (Smol) / (Ssfr) inside 0.75 r25 for each galaxy. Unless 
otherwise noted, we calculate J^mol using fixed oco and SsFR from 
Ha + 24^m. Quoted error bars on t^"°' report la scatter, uncertainties 
on the rank correlation arise from randomly repairing data. 



explore the effects of varying our approach to estimate 
SsFR and Smol- 

3.2.1. Choice of SFR Tracer 

Figure [2] and the lower part of Table [3] report the 
results of varying our approach to trace the SFR. We 
show SsFR estimated from only Ha, with a fixed, typi- 
cal Ana = 1 mag (top left) , along with results combining 
FUV, instead of Ha, with 24/im emission (top right). We 
also show the results of varying the approach to the IR 
cirrus. Our best-estimate Ssfr combines Ha or FUV 
with 24/im after correcting the 24/im emission for con- 
tamination by an IR cirrus following L12. We illustrate 
the impact of this correction by plotting results for two 
limiting cases of IR cirrus correction: no cirrus subtrac- 
tion (bottom left) and removing double our best cirrus 
estimate (bottom right), which we consider a maximum 
reasonable correction. Data density contours in Figure 
[2] weight each point equally and the large black points 
indicate the original binned results from Figure [1] 

The top left panel of Figure[2]and Table|3]show that the 
basic relationship between Ssfr and Smoi persists even 
when we derive Ssfr from Ha alone. The median t^I'^ 
and scatter using only Ha resemble what we find for our 
best estimate and the correlation between Ha and CO 
appears only moderately weaker than for the hybrid SFR 
tracer. It also appears moderately flatter than relations 
that incorporate IR emission as we underestimate ex- 
tinction in the central parts of galaxies fi j3.3l) . Inasmuch 
as Ha represents an unambiguous tracer of recent star 
formation, the top left panel in Figure demonstrates 
that subtle biases in the treatment of IR emission, e.g., 
24/im emission tracing the ISM rather than recent star 
formation, do not drive our results. 

The top right panel shows Esfr traced by FUV+24/im 
emission. The distribution agrees well with what we 
found using Ha+24/im, as do the median and scatter 



,_mol 
' dcp ■ 



The agreement of FUV+24/im and Ha with 
our best estimate Ha+24/.tm occur partially because we 
have designed our SFR tracers to yield self-consistent 
results (L12). However, that procedure considered only 
the overall normalization and did not require the detailed 
agreement we see comparing Figures [T] and [2l 

In the bottom row, we vary our approach to the in- 
frared cirrus. By default, we correct the 24/im map 
for infrared cirrus following L12. The bottom left panel 
shows the results of applying no cirrus subtraction, while 
in the bottom right panel we double our cirrus subtrac- 
tion. Turning off the cirrus subtraction yields median 



^mol 
'dcp 



10% shorter than our best estimate with notably 
lower scatter and our strongest observed correlation. The 
tighter correlation reflects the fact that the relationship 
between 24/tm and CO emission is the strongest in the 
data (see also Schruba et al. 2011). Esfr tracers that 
more heavily emphasize 24/tm exhibit the strongest cor- 
relation with Smol traced by CO. 
Doubhng the cirrus subtraction leads to an « 25% 



longer 



dcp 



larger scatter, and a mildly weaker correla- 



tion between Esfr and Smoi- This partially reflects un- 
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Fig. 2. — Star formation rate surface density, EgpR, estimated using different tracers as a function of molecular gas surface density, Snjoi. 
Light gray points show individual measurements, dark gray points indicate upper limits, contours show data density, and red points show 
the running median and standard deviation using each tracer. For comparison, black-and-white points show the bins from the top right 
panel of Figure [T] We show Ssfr estimated from: (top left) Ha assuming 1 magnitude of extinction, {top right) a combination of FUV 
and 2ifim emission (L12), (bottom left) Ha combined with 24fiTa emission with no IR cirrus correction applied to the 24/.im, (bottom right) 
Ha combined with 24/jm with double our IR cirrus correction applied to the 24/im. 



certainty in the cirrus calculation, which relies on model 
fits to observed data. It also reflects the deemphasis of 
24^m emission, which exhibits a very tight correspon- 
dence to CO, in favor of Ha, which still exhibits a good 
correspondence but with more scatter. The fraction of 
data that have upper only limits for Ssfr also increases, 
so that extending this analysis to lower surface density 
will require improved data and methodology. 

Thus we observe subtle variations in the relation be- 
tween EsFR and Smoi depending on the exact treatment 
of 24/im emission, including up to a ~ 50% variation in 
T^^°p across the full plausible range of cirrus treatments. 
However, our main results of a simple correspondence be- 
tween Sjiioi and SsFR hold even when we omit IR data 
from the analysis. Note that this conclusion relies on 
the assumption that Ha emission traces recent star for- 
mation. If a substantial fraction of Ha emission arises 



from sources other than recent star formation or if the 
mean free path of an ionizing photon regularly exceeds 
one of our kpc-sized resolution ele ments, then this gen- 
eral agreement may break down (j Rahman et al.l 120111 : 
ILiu et al.l[20Tll see discussion in L12). These more ex- 
otic situations aside, overall Figure [D and Table [3] show 
good qualitative and quantitative agreement among dif- 
ferent approaches to estimate Ssfr- We will find the 
same when fitting the data in Section 13.31 

Throughout the rest of the paper, we adopt Ha+24/im, 
corrected for the effects of a 24/im cirrus, as our single, 
best estimate of Ssfr- L12 justify this choice and we 
refer to that paper for more discussion. 

3.2.2. Choice of CO Lme 

HERACLES consists of maps of CO (2-1) emission, 
which we use to estimate the distribution of H2. CO (1- 
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Fig. 3. — Star formation rate surface density, SgFR, as a func- 
tion of molecular gas surface density, Emol, here estimated from 
literature CO(l-O) data (Section[2]|- Gray points show individual 
measurements, contours indicate data density, and red points show 
the running median and standard deviation using each tracer. For 
comparison, large black-and-white circles show the same bins from 
Figure m 

0) has been more commonly used to trace the distribution 
ofHaB. CO (1-0) maps exist for a s ubset of our targets 
(jHelfer et all [2001 IKuno et all [20071 ) . Though these do 
not have the same overall quality as the HERACLES 
maps, we use them to assess the impact of our choice of 
molecular gas tracer. Figure[3|plots Ssfr as a function of 
Smoi estimated fro m literature C O(l-O ) dat a. We allow 
repeat s, so that if iHelfer et al.l ()2003l) and iKuno et"al] 
(|2007l ) each mapped a galaxy we include each data set in 
Figure [H 

Overall, results for CO (2-1) and CO (1-0) agree fairly 
well. The CO (1-0) data tend to yield high er r^°i. This 
is exclusively a product of the 'Kuno et aLjdata. In the 
overlap of o ur sample and the Kuno et alJ data, the 
iKuno et"al1 data yield median r™°p 3.3 Gyr. Our 

data yield median T™°p « 2.3 Gyr for the same points. 
Howe ver, for the overlap with BIMA SONG (Helfc r et al.l 
l200l . the BIMA SONG CO (1-0) data give median 
2.1 Gyr. Over the same points, HERACLES 



,_mol 
'dop 



implies tJ^"' ~ 2.2 Gyr. The disagreement between our 
CO (2-1) data and CO (1-0) data thus appears no larger 
than the disagreement among published CO (1-0) data 
sets. 

3.2.3. CO-to-H2 Conversion Factor 

In Figures [1]- [3] we adopt a fixed CO-to-H2 conversion, 
aco • This assumption may be too coarse given the wide 
range of metallicities, dust-to-gas ratios, and central 
CO concentrations in our targets. Figure [4| shows Esfr 
and Smoi after the application of our "S = 100" con- 
version factor (Section I2.2p . which attempts to account 
for the presence of CO-poor envelopes of molecular 
clouds and central CO depressions. This conversion 



factor assumes that all CO emission originates from 
clouds with surface densities of 100 yi p, pc~^ with 
FDR structure like that described in iWolfire et al.l 
()2010( ) . The resulting dependence of aco on metallicity 
appro x imates the curren t cons e nsus (IWolfire et all 
20101: [(^er fc Mac Lo-^ [Mil ILerov et al.l l20Tir 



Feldmann et al.ll2012t iNaravanan et al.ll2012l) .~but note 



We emphasize that ISandstrom et all II2012I) demonstrate the 
ability o f CO (2-1) to robustly trace molecular gas in our sample 
(Section [22I1 



that this remains highly uncertain because of limited 
observational constraints. We also report results for the 
"E — 50" conversion factor in Table [3] This conversion 
factor makes the more aggressive assumption that a 
large amount of molecular emission emerges from weakly 
shielded parts of clouds, rendering aco very sensitive to 
the dust-to-gas ratio. 

Applying the "E = 100" aco to the data in Figure 
[1] yields Figure [4] The top rows, which show the bulk 
distribution of the data, appear qualitatively similar in 
the two plots, though Table [3| and close inspection of the 
plots do indicate that the normalization of the EsFR-Emoi 
relation changes between the two plots. The median 
weighting each line of sight equally rises from 2.2 Gyr to 
2.6 Gyr with the application of the variable conversion 
factor. 

The most dramatic contrast between Figures [1] and [4| 
appears in the bottom rows. Many of the low t^°' (high 

Esfr) data in Figure [T] arise from small galaxies with low 
dust-to-gas ratios. With the application of a variable 
conversion factor, our estimate of Smoi in these galax- 
ies moves to higher values while Ssfr remains constant. 
The result, visible in the bottom rows of Figure [U is that 
data from these low-mass, low metallicity galaxies now 
overlap the other points, forming a (more) continuous 
single SsFR-Sinoi trend. Table [3| reports that the scatter 
among galaxies drops from w 0.3 dex to ~ 0.2 dex with 
the application of the "S — 100" conversion factor, with 
the systematic difference in between high and low 
mass galaxies reduced from a factor of ^ 4 to a factor of 
- 2. 

Thus, as we will see in Section|4l application of a dust- 
to-gas ratio-dependent conversion factor CO-to-H2 con- 
version factor does affect the derived SspR-Smoi relation, 
with the sense of moving many points with low appar- 
ent r^p' into closer agreement with the distribution de- 
fined by large galaxies. This scenario of a rapidly vary- 
ing conversion factor and a weakly varying r^"' has been 
dis cussed in the c ontex t of the Small Magellanic Cloud 
by Bolatto et al. (2011) and in a theoretical context by 
[Krumholz e t al. (2011) • They interpret weak variations 
of r^p^ but strong variations of aco as evidence that the 
requisite preconditions for star formation more closely 
resemble those for H2 formation than those required for 
a high CO abundance. 

In 21 will show that the even more aggressive 
"E = 50" case may in fact explain most apparent galaxy- 
wide variations in r^gp , but note that while the "S = 50" 
conversion factor may offer an explanation for our ob- 
served trends, it also requires that a substantial fraction 
of the integrated CO emission from galaxies arise from 
relatively low Ay regions. 

The other portion of the "S = 100" conversion factor, 
the aco depressions in the central parts of galaxies, affect 
too few data to be prominent in Figured However, one 
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Fig. 4. — Star formation rate surface density, SgpR, as a function of molecular gas surface density, Smoh here estimated from CO (2-1) 
data combine d with our "S = 100" co nversion factor acOi which reflects variations in the dust-to-gas ratio and the central depressions in 
oco found bv lSandstrom et al.l II2012I ). Panels and annotations are as Figure [T] 

moi „^;„fc. o+ m^^cr-o+c j-jve best-fit relations for our ensemble of measurements 

and individual galaxies. We fit a relation with three pa- 
rameters: a normalization, a, power-law index, n, and in- 
trinsic, log-normally distributed scatter with RMS mag- 
nitude s. Then 



can see many individual low tJ^" points at moderate 



Smoi on close inspection. We return to this point in §43 

3.3. Power Law Index 

Studies of the star formation-gas connection in galax- 
ies have treated the relationship as a power law and fo- 
cused on the index of this power law. While this single 
parameter undoubtedly makes for easy shorthand, the 
fixation on this parameter obscures environmental fac- 
tors other than Egas. Recent observations offer good 
evidence that Ssfr relates to Shi and '^moi in funda- 
mentally differ ent ways (LOS, BOS iWong fc Blita [20021: 
iSchruba et 10120111) and tha t S SFR is even a multivalued 
function of S ^„i (iDaddreTal. 2010: Gcn zcl ct al. 2010; 



SsFR — a 



10 M© pc2 



(10) 



ISchruba et al1l2011l : lSaintonge et al.ll201lL l2012L 'and this 



paper) . 

Despite the shortcomings of this approach, we consider 
the best-fit index in our data as a useful, or at least ex- 
pected, point of comparison to previous studies. We de- 



with data intrinsically scattered by s. We derive the best- 
fit a, n, and s using a Monte Carlo approach based on the 
work of lBlanc et al.l ()2009[ ) . This resembles Hess diagram 
fitting used for optical color-magnitude diagrams. It in- 
cludes observational uncertainties, upper limits, and in- 
trinsic scatter in the relation. This approach also avoids 
important biases that can easily arise fitting scaling rela- 
tions to noisy, bivariate data. We illustrate these biases, 
which affect many commonl y adopted approaches, in the 
appendix (see also lBlanc et al...2009.) and note that they 
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Fig. 5. — Results of power-law fitting for (left) our combined data set and (right) individual galaxies. Left: Goodness of fit as a function 
of power law index, n, for fits to our entire en s emble of data where > 5 Mq pc~^. The fitting uses a modified version of the 

Monte Carlo technique outlined bv IBlanc et al.l I I2009I ). From top to bottom, we show results for fixed aco using our best Ifa+24/.im 
and FUV-|-24/im estimates; Ha-|-24/im with no cirrus removed or double the best-estimate cirrus removed; and If« with an assumed fixed 
i magnitude of extinction. We also show results for our "E = 100" and "E = 50" conversion factors. Vertical lines indicate where n reaches 
its minimum value, bold for the value in that panel, dotted for other panels. Each case shows a broad minimum in with a best-fit value 
near n 1 ± 0.2. Right: Best fit power-law index, n, for individual galaxies using Ha-f 24/jm and a fixed CO-to-H2 conversion factor. The 
distribution of best-fit n for individual galaxies peaks near n 1 with most galaxies in the range n 0.8-1.2 but outliers up to n > 1.5. 



TABLE 4 

Results of Monte Carlo Fitting to Equation [Tol 



Tracer 


log 




a 




n 


s 




(M© yr- 




kpc-2) 






(dex) 




-2.40 


± 


0.09 


0.95 


±0.15 


0.30 ±0.05 


... no cirrus removed 


-2.35 


± 


0.08 


0.95 


±0.13 


0.225 ± 0.05 


... double cirrus removed 


-2.50 


± 


0.10 


1.00 


±0.17 


0.375 ± 0.05 


FUV-f24/im 


-2.40 


± 


0.09 


0.95 


±0.16 


0.30 ±0.05 


Ha with Ajja = 1 mag 


-2.35 


± 


0.09 


0.80 


±0.14 


0.275 ± 0.06 


S = 100 Qco and Ho+24/im 


-2.45 


± 


0.09 


1.0 


±0.15 


0.30 ±0.05 


S = 50 aco and Ha4-24^m 


-2.35 


± 


0.09 


0.90 


±0.15 


0.30 ±0.05 



Note. — Results of Monte Carlo fitting to Equation 1101 for different combi- 
nations of EgFR and En^oi tracers. Column (1) reports the tracer used; in the 
top five rows we vary the Ssfr tracer while adopting fixed aco- The last two 
rows adopt our best Esfr estimate, Ha+24/im, and vary the adopted conversion 
factor. Columns (2)— (4) report the best-fit coefficient at Snjoi = 10 Mq pc~^, 
the power law index, and the intrinsic scatter. We quote uncertainties from the 
Monte Carlo simulations described in the appendix. 



can easily shift the derived index by a few tenths for 
reahstic data distributions. 

Following IBlanc et all (pOOl ) we grid our data, deriv- 
ing a two-dimensional image of data density in regularly- 
spaced cells in l ogi „ SsFR-logio ^^moi space. Unlike 
IBlanc et al.l (|2009t ) we work in logarithmic space. This 
gives us a better ability to resolve the distribution of 
our data, but forces a coarser approach to upper limits. 
We treat upper limits by essentially creating an "upper 
limit row" along the Ssfr axis. In detail, we adopt the 
following approach: 



1. We exclude all data with Smoi < 5 M© pc ^. This 
gives us a data set with a well-defined x-axis. 

2. We generate Monte Carlo data sets for a wide range 
of a, n, and s in the following way. We take our 
observed Smoi to represent the true physical distri- 
bution. We draw 100,000 data points from this dis- 
tribution (allowing repeats) for each combination 
of a, n, and s. We derive Ssfr for each of these 
points. We then apply the expected uncertainty 
to Smoi (the statistical uncertainty from HERA- 
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CLES) and Ssfr (0.15 dex). We grid these data in 
logiQ EsFR-logio Snioi space, using cells 0.125 dex 
wide in both dimensions. We treat this grid as the 
expected probability distribution function for those 
underlying parameters a, n, and s. 

3. We grid our observed data in logj^g SgpR-logj^o ^^moi 
using the same grid on which we derived probability 
distribution functions. We create different grids for 
each set of Ssfr estimates. 

4. We compare our gridded data to the Monte Carlo 
realization for each a, n, s combination and cal- 
culate a goodness-of-fit estimate, which we here 
loosely refer to as x^- After re- normalizing the 
Monte Carlo grid to have the same amount of data 
as the observed grid, we calculate: 



E 



(^obs - ^modcl)' 



(11) 



model 



where the sum runs across all cells, i, iVobs refers 
to the observed number of data in the grid cell, 
and A^modci refers to the expected number of data 
in that cell given a, n, and s and our observational 
uncertainties. The goodness-of-fit statistic is thus 
analogous to calculated for the case of Poisson 
noise in each cell. Points with only upper limits 
on SsFR (where Ssfr < 10"'^ Mq yr^^ kpc~^) are 
included in the calculation. These have an associ- 
ated Smoi value but all upper limits are treated as 
having the same Ssfr- 

We apply this method to our ensemble of data, repeat- 
ing the exercise for each SFR tracer discussed in Section 

and for our fixed, "S = 100" , and "S = 50" conver- 
sion factors. We also fit each galaxy on its owrF^. Figure 
[S] and Table |3] report our fits to the combined data set. 
Figure [5] plots the approximate reduced a function 
of power law index, marginalizing over a and s. We ob- 
serve clear minima in the range n = 0.5-1.5 for all SFR 
tracers. The appendix presents a Monte Carlo treat- 
ment that considers a number of effects: robustness to 
removal of individual data or galaxies, statistical noise, 
calibration (gain) uncertainties for each data set, and ex- 
act choice of fitting methodology. We quote uncertainties 
derived from this Monte Carlo treatment in Table 01 

The fits in Table 2] suggest a power law with n w 

1 ± 0.15, intrinsic scatter of a factor of ~ 2 (0.3 dex), 
and T™°i « 2.5 Gyr at E^oi = 10 Mq pc'^. The 
slope remains consistent wit h a linear relation between 
H2 and star formation (BOS, IBigiel et al.ll2011l) or with 
the weakly super-l inear slope of Genzel et al.l ()201ft) or 
iDaddi et al.l ()2010D . though note that our 1-kpc scale 
does not precisely match their observations. The mild 
difference between the best-fit coefficient, a, and the me- 
dian T^p' reported in Table [3] reflect the inadequacy of 
the power law to capture the full distribution of the data. 

The choice of SFR tracer affects the fit, but offers more 
of a refinement than a qualitative shift in these conclu- 

Due to the lower density of points, we use 0.2-dex cell sizes and 
require only 10, 000 points to populate the theoretical distribution. 



sions. The sense of the shifts resemble those seen in Sec- 
tion l3.2l Replacing FUV for Ha as the unobscured tracer 
has minimal effect. Using only Ha to estimate SFR yields 
a slightly shallower slope. Based on the observed Ha-to- 
IR ratio, exti nction inc reases with increasing Ssfr (see 
plots in .Prescott et al.l [2007. and L12). By assuming a 
fixed Ayia we would expect to underestimate Ssfr at 
the high end and overestimate it at low end, somewhat 
"tilting" the relationship to shallower slope. If we do not 
remove any cirrus from the data, working only with the 
measured 24/xm emission, the scatter in the relation di- 
minishes to less that 0.2 dex. This underscores the point 
that it is the tight observed correlation between CO and 
24/xm emission that drives much o f the recent work on 
this topic (see mor e discussion in iRahman et all 120111 : 
ISchruba et al. 2011: ILiu et"aIll201lD so that SFR tracers 
that emphasize 24/im data tend to yield the tightest re- 
lations. Conversely, increasing the cirrus removed leads 
to a somewhat longer overall 7"^^^°' with larger intrinsic 
scatter. Adjusting the conversion factor exerts only a 
mild impact on the fit because largest corrections apply 
to small galaxies and often to low apparent Smoi regions. 
These significant variations to a small subset of the data 
do not drive substantial variations in the fit. 

The simple nearly linear scaling given by our fits could 
result from the superposition of a varied set of distinc t 
relations for individual galaxies (seel Schruba et al.l201lD . 
In the right hand panel of Figure [5] we show that indeed 
the best-fit index for individual galaxies exhibits signif- 
icant scatter. We find a median 1.05 but best-fit values 
span ~ 0.5-1.5. We report best-fit indices for individual 
galaxies in the appendix and stress two general conclu- 
sions here. First, we see variation in index from galaxy 
to galaxy, but the 67% range is still 0.8-1.25, consistent 
with the idea that to first order the molecular gas supply 
regulates the star formation distribution and in sharp 
con trast to the steep i ndices relating Ssfr, to atomic 
gas IBigiel et al.l (|20Cl8f) : ISchruba et all (|2011D . Second, 
the fact that these galaxy-to-galaxy variations wash out 
into Figure [1] implies that while t^I'^ may correlate with 
Smoi within an individual galaxy, simply knowing Smoi 
at 1 kpc resolution with no other knowledge of local con- 
ditions or host galaxy does not allow one to predict t^^'^ 

better than simply adopting a median t^"'. That is, in 
the absence of knowledge of other conditions, Smoi is not 
a good predictor of the molecular gas depletion time. 

3.4. Comparison to Literature Data 

Many studies have assessed the relationship between 
gas and star formation in nearby galaxies (Section [T]) . 
Figure [S] and Table O compare our measurements to 
a compilation of these studies (see Section 12.5! and 
IBigiel et al.ll201 J ). We adjust each set of measurements 
to share our adopted CO-to-H2 conversion factor and 
stellar initial mass function. 

Table [5] give s r^"' by study. These span from 0.4 for 
the lKennicuttI (Il998bf) starbursts to 4.4 Gyr f or the study 
of NGC 6946 bv lCrosthwaite fc Turned ((200I . Consider- 
ing all measurements equally, the median literature 
is 2.7 Gyr, whi ch drops to 2.1 Gyr if we exclude the 
large data set of IRahman et al.l (|2012 f). which otherwise 
dominates the statistics. Treating each study as a single 
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independent measurement, the median is r^"' = 2.0 Gyr. 
These are in good agreement with the estimates of this 



study (Table El) 



or 

'dcp 



2.2 Gyr weighting ah measure- 



FlG. 6. — SsFR versus Ejnol for our results and literature measurements. The type of literature measurement varies with panel: the top 
left shows galaxy averages, the top right shows radial profiles, the lower left shows sparse, coarse (~ 1') sampling of galaxy disks, and the 
lower right shows fine sampling (< 15"). In all four panels, red points show a running median and standard deviation (error bars) for our 
kpc-scale measurements (Figure[IJ using Ha+24/^m and CO_(2-l). Orange circles in the top left panel show integra ted measurements for 
our data. Literature measurements are describe d in Se ction references: K98 — 'Kennicutt? fl998a), CIO + lit. — ICalzetti et al.l (20W[) 
combined with litera ture C O, S12 — Saintong e et al.l 1201 2) with offset to bring SFR estimates i nto agreement (S ection I2.2[l. W02 — 
IWong fc BTItg 1 I2OOI). SOT — ISchuster et al., t2007l). C07^Crosthwaite &: Turner. (2007.1. M02 — .Murgia et al.l 12003'), L05 — ILerov et aU 
H2005I ). KoT^ IKennicutt et al.l H2007I ). B09 — IBIanc et al.l H2009I 1. Rll — IRahman et al.l l(20Tll ). R12 — IRahman et al.l I^OTW . 

60". The bottom right panel shows studies that obtain 
high-angular-resolution sampling of each target. In each 
panel we plot the running median and standard deviation 
for our data, binned by Smoi, as red points. 

The final three panels of Figure |6] demonstrate excel- 
lent agreement between our data and previous studies 
that resolve the disks of galaxies (see also Bigicl et al] 
l2011f ). This agreement may not be surprising given 
that our study shares targets with many of these lit- 
erature studies, which also heavily overlap one another. 
Nonetheless, we show here that repeated measurements 
of the distribution of Ssfr and Smoi in the nearest star- 
forming spiral galaxies mostly cover the same part of 
parameter space regardless of exact methodology. Un- 
certainty in interpretation and fitting techniques have 



ments equally. The scatter among individual literature 
data is ~ 0.36 dex and from study to study the scatter 
is « 0.23 dex. 

Figure [S] shows a more detailed comparison between 
our measurements and individual literature data. We 
separate the literature studies according to the scale sam- 
pled. The top left panel shows measurements where one 
point corresponds to one galaxy. The top right panel 
shows data from studies that measure azimuthally aver- 
aged Smoi and SgpR in a series of concentric tilted rings. 
The bottom left panel shows data for individual point- 
ings with comparatively poor angular resolution, « 40- 
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TABLE 5 

Comparison to Literature 



Galaxy 


Median r™°' 

dcp 


Scatter 




[Gyr] 


[dex] 


This study ... 






... weighting by measurement 


2.2 


0.28 


... weighting by galaxy 


1.3 


0.32 


Median of literature . . . 






... weighting by measurement 


2.7=* 


0.36 


... weighting by study 


2.0 


0.23 


Uislt Averages 


Saintonge et al. (2012) 


0.7 


0.37 


... offset SFR^ 


1.1 


0.37 


Calzetti et al. (2010") + literature CO'' 


1.5 


0.38 


Kennicutt (1998^7".. 






... disks 


1.1 


0.46 


... starbursts 


0.4 


0.39 


Radial Profiles 


Schuster et al. (2007) 


2.0 


0.12^* 


Crosthwaite & Turner (2007) 


4.4 


o.oe'^ 


Wong & BUtz (2002)" 


2.0 


0.36 


Low Resolution 


Lcrov et al. (2005) 


2.1 


0.33 


Murgia et al. (2002) 


2.8 


0.41 


High Resolution 


Rahman et al. (2012) 


2.9 


0.37 


Rahman et al. (2011) 


1.6 


0.15^^ 


mane et aL (2009) 


3.2 


0.64^^ 


Kennicutt et al. (2007) 


2.2 


0.37^^ 



Note. — Average molecular gas depletion time, in Gyr, for 
matched assumptions — ^CO = 2 X 10 cm~^ (K km s~^)~^, a 
Kroupa IMF, and including helium in the gas estimate. The left 
column gives the study, with the list broken down by sampling 
approach and the right column reports the median molecular gas 
depletion time in that study. Error bars report the Icr scatter, in 
dex, for each study. 

^ Dominated by IRahman et al.l II2012I ) . Without I Rahman et al.l 
(12012) median is 2.1 Gyr. 

^ SFR estim ate offset to matc h our estimates (Section 12.21 1. 
" CO from Youne et al.l ||T995|) and Helfcr ct al. ( 2003). 

Study considered a singe galaxy. Others combine multiple galax- 
ies. 

clouded this basic agreement in where the data lie. 
Given our basic approach to physical parameter estima- 
tion, there appears to be overall agreement for a typical 



,_mol 
'dcp 



1-3 Gyr in local disk galaxies. 



The first panel of the Figure IH] looks qualitatively dif- 
ferent from the other three. This panel shows galaxy- 
integrated measurements, so that one point is one galaxy. 
These tend to scatter from overlapping our data up to 



lower 



,_mol 
'dcp 



at low Emoi and comparatively high Ssfr- 



This same effect appears in Table [5] as low values of r|^p 
for studies that focus on measurements of whole galaxies. 
Our synthesis of literature CO and SFR measu r ements 
yields T^°i = 1.3 Gyr while for the iKennicut^ (|1998al ) 

disk galaxies the median r^^"' = 1.1 Gyr. Treating our 
own sample as a set of integrated measurements we find 
a similar value, r^"' — 1.3 Gyr. 

This disk-integrated r^^"' is significantly shorter than 
the that we measure treating each point equally. We 
noted this effect in Section [3Tl It arises because weight- 
ing each galaxy equally emphasizes small, low-mass, low 
SFR galaxies relative to large, massive galaxies. These 
low mass galaxies have less physical area than large disks, 
so that they do not affect the ensemble of measurements 



much. However these small galaxies do exhibit short ap- 
parent rj^p' and when given equal weight they drive the 
median r^"' down by a factor of ^ 2. We explore this 
and other systematic variations in t^I'^ in the second part 
of this paper. 

4. SYSTEMATIC SECOND-ORDER VARIATIONS IN rMpL. 
GLOBAL CORRELATIONS, EFFICIENT GALAXY 
CENTERS, AND CORRELATED SCATTER 

In Section [3] we demonstrate that our ensemble of data 
can be described to first order by a roughly linear relation 
between Ssfr and T^moi with a slope corresponding to a 
typical depletion time t^'^ 2.2 Gyr with a factor of 
two scatter from line of sight to line of sight. However, 
we also show that the apparent uniformity of t^^"' results 
at least partially from the emphasis that our approach 
places on the disks of large, star-forming galaxies. These 
contribute most of the area in our sample. When we 
apply weightings that emphasize small galaxies or the 
inner parts of galaxies, we observe departures from this 
simple picture. 

In this section, we explore these variations. We ex- 
amine the dependence of disk-average r^^"' on integrated 

galaxy properties and the dependence of r^"' on 

local physical conditions (i i4.2p . We highlight apparent 
variations in t^"' as a function of galaxy mass, metallic- 
ity, and dust-to-gas ratio and we discuss the CO-to-H2 
conversion factor as a potential cause. In Section we 
contrast the central regions of our targets with the disks 
and show evidence for a shift to more efficient star forma- 
tion in galaxy centers, perhaps indicative of a transition 
between "disk" and "starburst" modes of star formation. 
Finally we examine the scale-dependence of scatter in 

mol 
Mcp 



to show that undiagnosed systematic varia- 



tions in Tjlj^p persist in our data, reflecting either large- 
scale synchronization of star formation or real correlated 
efficiency variations spanning the disks of our targets. 



4.1. 



Galaxy-to- Galaxy Variations in r^""' 



Our comparison with the literature reveals the same 
effect that we saw in Section 13.11 when we weight by 
host galaxy instead of by individual kpc-scale measure- 
ment, clear variations in t^^"' become visible. In par- 
ticular, a large population of short r^"' data emerges 
in the bottom left panel of Figure [1] and the top left 
panel of Figure [51 In Figure [7] we investigate the phys- 
ical origin of these differences. Adopting a fixed aco-, 
we plot T^'^ for each galaxy in our sample as a func- 
tion of a series of galaxy properties: stellar mass, disk- 
average total (H I -I- H2) gas surface density, metallicity, 
dust-to-gas ratio, rotation velocity, average orbital time 
(weighted by Smoi), morphological type, and specific star 
formation rate (sSFR=5'i^i?/M,). r^°' correlates with 
many of these quantities with the overall sense that low 
mass, low metallicity, late-type, high sSFR galaxies ex- 
hibit shorter apparent molecular gas depletion times than 
massive, metal-rich, earlier type disk galaxies. 

Such trends have been observed in various ways be- 
fore. In Figure [3 gray points show the appr oximate r^"' 
as a function of morphology as inferred by lYoung et al.l 
(|1996i) from comparing the FCRAO Extragalactic CO 
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Fig. 7. — Integrated (Sjhoi/Ssfr) within 0.75 r25 and standard deviation (error bars) in T^°p over individual lines of sight plotted 

against integrated galaxy properties: stellar mass, average gas mass surface density, metallicity, dust-to-gas mass ratio, rotation velocity 
at r25, orbital pe riod at 0.4 r25 , morpholo gy, and spe cific star formation rate ((Sgppj^) / (2*)). Literature data appear in gr ay (Y96, Kll, 
and Sf2 refer to I Young et al.| [T996: Krum holz et al.| i201f : Saintongc ct al. 2012) and fits or predictions appear as curves: IWolfire et al.l 
II2010I ) conversion factors for E = 50 and 100 IsAq pc~^, the rotation- velocity threshold associated with the suppression of dust lanes by 
IDalcanton et all f2004, D04), a fixed efficiency per orbital period, and the t^°' vs SSFR fit by ISaintonge et al.l i'li) I ll Sll). Low mass, 
metal poor, late type galaxies show significantly lower 1"^°' (high SFR-to-CO ratios) than their massive, early type counterparts, but these 

trends persist with smaitcr dynamic range even among massive gaiaxies (we mark galaxies with Af, > lO"*^" Mq using magenta dots in all 
panels). 

found a clear relation between t^I'^ and Af* and iden- 



Survey with Ha (squares) and IR (diamonds) data. They 
found later-type galaxies to exhibit mildly enhanced star 
formation efficiencies. We also plot t^'^ derived by 



iKrumholz et al.l ()201lD combining literature SFR and 
CO data (including HERACLES measurements), who 
showed a strong trend of S FR/CQ decreasing w ith in- 
creasing metallicity (see also lSchruba et ani2012f) . Both 
literature trends agree well with our own data and here 
we present the first direct comparison to the dust-to- 
gas ratio, the quantity thought to link metallicity to 
both CO emission and the star formatio n process (e.g., 
IWolfire et aLllMot iKrmnholz et al.l[20Tll). 
ISaintonge et all ()201lD and ISaintonge et al.l ()2012l ) 



tified sSFR as the global galaxy property most directly 
linked to Tj"°'. We plot their data, corrected to bring our 
SFR estimates into agreement in the top row and bot- 
tom right panel. In the bottom right panel, we plot their 
fit of T^^p , corrected for diflFerent treatments of helium. 
We al so indicate the r otation velocity threshold below 
which lDalcanton et al.l (|20M) noted the disappearance of 
prominent dust lanes in edge-on disk galaxies. The exis- 
tence of such dust lanes should indicate the presence of a 
dusty, dense ISM of the type that mig ht host molecular 
material and star formation. Though IDalcanton et al] 
(|2004) make no clear prediction for t^"' we note that 
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Fig. 8. — Same as Figure [7] but after application of our variable "E = 100" CO-to-H2 conversion factor. Most of the galaxy-scale 
variations in "rj^"' visible in Figure [7] can largely be explained by this dust-to-gas ratio- dependent conversion factor, though significant 
galaxy-to-galaxy scatter remains. A dashed horizontal line in each panel shows the median t^"°' across our sample and we reproduce several 
comparison lines from Figure [7] 



there is a generally lower t^'^ below the threshold than 
above; at the very least, the lack of abundant dust shield- 
ing may d rive both observation s. This contrasts with the 
finding of iWatson et al.l ()2012|) who do not find such a 
trend with circular velocity. 

Table |6] quantifies the relationship between t"1'^ and 
these other quantities using the Spearman rank correla- 
tion coefficient. t^°' correlates with stellar mass, metal- 
licity, dust-to-gas ratio, and rotation velocity at > 3a 
significance. Weaker (~ 2a) correlations link t^°' to 
gas surface density, orbital time, morphology, and spe- 
cific star formation rate. The trends with stellar mass 
and rotation velocity appear continuous, without clear 
thresholds. Indeed they still emerge, though at lower 
significance, even if we consider only high mass galaxies. 

The sense of the correlation with gas surface density 
is that we find the shortest depletion times for galaxies 



below (Shi+H2) ^ 10 M0 pc^^. These systems still yield 
the lowest t^'^ regardless of our treatment of aco (see 
below). In our sample, these 1ow-(I]hi+h2) hold a more 
of their gas in the atomic phase than systems with high 
(SHI+H2)- Therefore, the trend we see amounts to longer 
depletion times in gas-rich, molecule dominated systems. 
However, note that our high (SHI-1-H2) systems are gas- 
rich spirals. Our sample does not include truly extreme 
systems like the local U/LIRGs. Many of these systems 
have high gas surface densities but very low r^"' (e.g., 
iSanders &: Mirabeilll996[ ). so that if they were added to 
Figure [3 they would not extend the trend that we see. 
We do see similar efficient starbursts within our targets 
( §4.4p . but the effect remains mostly confined to the cen- 
tral parts of galaxies and does not propagate to these 
galaxy-scale measur ements. This behavior m ay be some- 
what evident in the iSaintonge et "all (|2012l) sample, but 
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TABLE 6 

Rank Correlation of t",""' with Galaxy-Average Properties 



Quantity Fixed «co ^ = 100 ctCO S = 50 «co 



Stellar mass 




+0. 


,72 ± 


0, 


,18 


-1-0.47 ±0.19 


-0.07 + 0.19 


... M. > 10^° 


M© 


n 

— u. 


'i7 -\- 
O 1 ^ 


n 
u. 


9n 
zu 






(ShI+H2} 




+0. 


,40 ± 


0, 


19 


+0.64 lb 0.18 


+0.49 ±0.20 


... M. > IQio 


Me 


-un 

-ru. 


19 ± 


Q 


20 


in fin -1- n 99 


in fin n 9^^ 


Metallicity 




+0, 


,73 ± 


0, 


,21 


4-0.44 ±0.19 


-0.14 + 0.19 


... M. > IQlO 


Mq 




59 ± 


0. 


20 


yj.±\j _i_ yj.^'j 


_n 91 + n 99 


Dust-to-gas ratio 


+0, 


,81 ± 


0, 


19 


+0.28 ±0.19 


-0.49 + 0.19 


... M. > IQlO 


Mq 


+0, 


,57 ± 


0, 


,20 


-0.24 ±0.22 


-0.72 + 0.23 


Rotation velocity 


+0, 


,74 ± 


0, 


,20 


+0.47 ±0.20 


-0.11 + 0.23 


... M. > 10^° 


M© 


-0, 


,40 ± 


0, 


,24 


+0.02 + 0.026 


-0.24 + 0.23 


Orbital time 




-0, 


,45 ± 


0. 


,21 


-0.36 + 0.19 


-0.12 + 0.20 


... M. > IQl" 


Mq 


-0, 


,32 ± 


0, 


,24 


-0.21 + 0.23 


-0.03 + 0.22 


Morphology 




-0. 


,49 ± 


0. 


19 


-0.30 + 0.19 


-0.03 + 0.20 


... M. > 10"' 


M0 


-0. 


,16 ± 


0. 


,23 


+0.03 + 0.23 


+0.16 + 0.22 


SFR/M, 




-0. 


,38 ± 


0, 


19 


-0.06 + 0.19 


0.32 + 0.19 


... M, > 10"' 


Mq 


-0, 


,27 ± 


0, 


,22 


+0.16 + 0.23 


0.34 + 0.25 



Note. — Rank correlation coefficient relating the molecular gas 
depletion time, t^°^ to galaxy-average properties. Quoted uncer- 
tainties report the Icr scatter of the correlation coefficient about 
obtained by randomly repairing the data. The columns give re- 
sults for different assumptions about the CO-to-H2 conversion fac- 
tor. We report results for all galaxies and only high mass galax- 
ies, M« > 10^" Mq. The average molecular gas depletion time 
is strongly covariant with galaxy average properties. This covari- 
ance can be reduced but not removed by the application of a D/G- 
dependent CO-to-H2 conversion factor. 

not our data. 

The relation between sSFR and r"^°' does not appear 
as strong as other trends in our sample, contrary to the 
finding of iSaintonge et aTl ()2011[ ). This may simply re- 
flect our sample's lack of massive early type galaxies with 
low sS FR or U/LIRGs with h igh sSFR. Where we do 
overlap ISaintonge et al.l ()2012|) . the agreement between 
the samples appears very good after our adjustment for 
different approaches to SFR estimation {^^. The most 
substantive obser vational disagreement between our re- 
sults and those of ISaintonge et ahl (|2011[ ) is that we find 
T^p to correlate with metallicity and dust-to-gas ratio, 
even for relatively high mass galaxies, while their tests 
revealed no such correlations. 

Dust-to-Gas Ratio and Conversion Factor: We em- 
phasize that the correlations in Figure [7] relate apparent 
'''Sep ^'^ integrated galaxy properties. The second row 
of Figure [7] shows strong trends in t"^!,"' as a function of 
metallicity and dust-to-gas ratio and these trends raise a 
fundamental issue regarding the interpretation of Figure 
[71 There is good evidence that the CO-to-H2 conversion 
factor increases with decr easing metallicity and dust-to- 
gas ratio (see refere nces in lIsraellFl997l : iLerov et"al1l2011l : 
iBolatto et al.l [20Tll ). Metallicity and dust-to-gas ratio 
vary with the stellar and dynamical mass of a galaxy, so 
that the correlations in Figure [7] and Table [6] may reflect 
either true variation in the efficiency with which molecu- 



lar gas forms stars, variations in the CO-to-H2 conversion 
factor, or a mixture of the two. 

The blue and red curves in the second row show the 
dependence of aco on metallicity or dust-to-gas ratio 
for the "S = 100" (blue) and "S = 50" (red) models 
described in Section [5] assuming a fixed molecular gas 
depletion time. The third and fourth columns in Ta- 
ble E] report the correlations between r^""' with global 
properties after application of these conversion factorj^. 
The red curve, representing the more extreme assump- 
tion that CO emission arises from low surface density, 
E = 50 M0 pc~^, clouds fits the dust-to-gas ratio trend 
well. Neither trend is steep enough to fit the metallicity 
data (left panel, second row), but given the large uncer- 
tain ties in the metallicities (often several 0.1s of a dex, 
e.g.. lMoustakas et al.|[20Tol ). the "E = 50" curve is not a 
bad match. 

Comparison of these two curves and the second row of 
Figure [7] highlights several technical points also found in 
the recent literature. First, the conversion factor needed 
to impose a fixed t^"' depends more sharply on metal- 
licity than the predictions from many theoretical mod- 
els a nd some observa t ional determinations (see plots in 
Kru mholz et al.ll2011l : iSchruba et alll2012l : iGenzel et afl 
l2012f ). Second the left and right panels of the second 
row do not perfectly agree. The dust-to-gas ratio de- 
rived from Spitzer observations of low metallicity regions 
appears somewhat lower than one would predict from a 
linear scaling of their metallicity (jMuhoz-Mateos et alj 
|2009( ). Because both the observed and predicted trends 
flatten at high D/G and high metallicities, the zero 
points (solar values) used to compute the theoretical 
curves might be adjusted to yield better agreement be- 
tween metallicity and D/G or to better match one the- 
oretical curve or the other to the data. That is, given 
the uncertainties in the absolute measurements of both 
metallicity and D/G, one can "slide" the relative posi- 
tions of the theoretical curves and the D/G and metal- 
licity determinations left and right. As a result, the dif- 
ferences between the "E = 50" and "S = 100" curves 
should not be overemphasized. 

These details aside, the key questions are: to what de- 
gree can conversion factor variations explain the global 
trends seen in Figure [7] and are the adopted conversion 
factors reasonable? In Figure [8] we repeat Figure [7] but 
adopt our "E = 100" conversion factor rather than fixed 
aco- Many correlations of r^"' with local properties ap- 
pear weaker with the "E — 100" aco than with a fixed 
conversion factor. The more extreme "E = 50" case re- 
moves almost all variation with stellar mass and metallic- 
ity, but reverses the trend with dust-to-gas ratio. Thus, 
the "E = 50" may represent the simplest single explana- 
tion for the bulk of the variation in Figure [T] "E = 100" 
corresponds to a conservative aco estimate that refiects 
present conventional wisdom. It lessens the strength of 
many correlations but (marginally) significant variations 
in T^p' with other properties still exist when using this 
conversion factor. 

Note that the conversion factors arc calculated from the dust- 
to-gas ratio profiles and then applied to the CO maps to derive 
the integrated H2 masses and that these conversion factors include 
central aco depressions, so that the mapping of conversion factor 
to the average dust-to-gas ratio is not perfect. 
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4.1.1. Conversion Factors and ISM Structure 

Both of our conversion factors assume a simple, uni- 
versal structure for the molecular ISM beneath our kpc 
resolution. In actuality, the typical surface density of 
GMCs, the universality of this value, and the balance 
of GMCs and more diffuse molecular material remain 
poorly constrained by observations. Much previous work 
places typical surface densities of Milky Way GM Cs at ^ 
150Mq pc~^ (Solomon_et_al. 1987; Ro man-Duval et al.l 
|201C1( ). High resolution observations of the inner parts 
of nearby spirals (NGC 5194, NGC 6946) y ield similar 
surface densities ([Donovan Meyer et al.ll2012l Hughes et 
al., in review; Colombo et al., submitted). Other stud- 
ies find lower surface densities, even in th e same sys- 
tems: e.g., studies of the Milk y Way ring bylHever et al.l 
dSooi) , the disk of NGC 6946 fe eboUedo et al.ll2012Ohe 
LMC (|Fukui et alllW)! IHughes et al.ll2ninl:IWong et al ' 



2 increase in r™°' moving from low 



2011 ). or a sample of Local Group galaxies ([Bolatto et al 



all find typical CMC surface densities E ^ 50- 
100 M0 pc~^. Some of these differences may be at- 
tributed to methodology, but Hughes et al. (in review) 
decisively demonstrate that the structure of the molecu- 
lar ISM does vary substantially with environment via a 
carefully matched comparison M51, M33, and the LMC. 
They show substantial differences in the probability dis- 
tribution function of CO emission at ^ 50 pc resolution 
both within and among these galaxies. 

Given these uncertainties, the best way to read the 
"S = 50" conversion factor is a case where aco depends 
sharply on the dust-to-gas ratio because most CO emis- 
sion to come from low-extinction lines of sight. Ay ~ 1- 
2 at Milky Way dust-to-gas ratios. In addition to the 
mixed evidence on surface densiti es of whole populations , 
maps of very loc al clouds (as in lHeiderman et all 120101 : 
iLada et aLll2010D do tend to find significant mass in low- 
Ay components. Furthermore, the c ontribution of dif- 
fuse sight lines (e.g., those studied bv 'Liszt et al.l[20To[ ) 
to the integrated CO emission from a large part of a 
galaxy remains poorly known, but may easily be several 
10s of percent in the Milky Way (H. Liszt, private com- 
munication) . 

A detailed investigation of CMC structure is beyond 
the scope of this paper. The key points are the follow- 
ing: (1) the basic structure of the molecular ISM remains 
uncertain at a level that significantly affects our interpre- 
tation of integrated CO emission, (2) the structure of the 
molecular ISM does vary substantially with environment, 
and (3) observations do appear to admit the possibility 
of substantial CO emission from low- Ay lines of sight or 
unbound clouds. Given these caveats, we must view our 
adopted conversion factors as more schematic than ex- 
act, so that subtle differences between the "E = 50" and 
"E = 100" case should not be over-interpreted. We can 
say from Table |6] and Figure [7] and [8] that the conversion 
factor can explain many of the largest systematic vari- 
ations in r|^p across our data, but probably not all of 
them. High mass, high surface density galaxies probably 
are less efficient at forming stars from their molecular 
reservoir (i.e., high t^I'^) than low mass galaxies. 

4.1.2. Low Efficiency (High t^"^) in High Mass Galaxies? 

Our best guess is that weak correlation exists relating 
■^dep to galaxy mass and average surface density, reflect- 



ing a factor of ^ 2 increase in t^'^ 

(M, < 10^" Mq) to high {AU > Mq) mass galaxies. 
Several natural explanations exist for such trends and a 
combination of these are almost certainly at play. Sup- 
pression of star formation in molecular gas by dynamical 
effects appears evident from high resolution observations 
of M51 (Meidt et al. submitte d) and in some ba rred spi- 
ral galaxies (see discussion in lJogee et al.ll2005h . These 
large-scale dynamical effects may suppress the ability 
of molecular gas to collapse into gravitationally bound, 
star-forming clouds. 

This represents a subset of a more general ef- 
fect: in regions of high pressure and high gas sur- 
face density, the ISM bec omes increasingly molecular 
(jBhtz fc Rosolowskvl [2001 LOS). A diffuse molecular 
ISM, meaning unbound material and the low density out- 
skirts of bound clouds, may represent a reservoir of non- 
star-forming molecular gas (though note that even the 
bulk of bound, molecular material does not d irectly par- 
ticipa te in star-formation; see references in iLada et al.l 
[20T1 . In low mass, low molecular fraction galaxies such 
gas, if present, will often be in the atomic phase. Higher 
mass galaxies will have both a higher molecular fraction 
and stronger dynamical effects such as shear and stream- 
ing moti ons working at suppress ing the creation of bound 
gas. See iSaintonge et al.l ([2012| ) for discussion of analo- 
gous effects. 

This schematic picture leads to several predictions and 
is clouded by several complications. Both a prediction 
and a complication is that in this sketch, the physical 
state of molecular gas varies systematically as a func- 
tion of galaxy mass. This should manifest via several 
direction observables such as the CO-to-HCN or CO-to- 
HCO-f ratio, which will trace the fraction of dense molec- 
ular gas, the ratio of ^■^CO-to--^^CO, tracing the optical 
depth of the gas, and the ratio of CO rotational tran- 
sitions, which trace a complex combination of density, 
temperature, and deviation from local thermodynamic 
equilibrium. An immediate complication from such vari- 
ations is that changes in internal conditions propagate to 
aco variations that are distinct from the dust-shielding 
effects accounted i n our adopted conversio n factors (e.g., 
see discussion in iMalonev fc Blackl Il988f ). A second, 
weak prediction would be a general increase in t^"' with 
increasing molecular fraction, at least up to a certain ex- 
tent — starbursts exhibit both low r^"' and high molec- 
ular fractions. Systematic trends in t^'^ with molecular 
fraction were not immediately clear in LOS; we search 
for them in our present data set in the next section. A 
similar mixed picture arises from dynamical effects: one 
might expect to see high molecular fractions but low ef- 
ficiencies in regions where shear or streami ng motions 
suppr ess bound cloud formation. However, iFovle et al.l 
([2OIOI ) found little or no evidence for t^'^ variations be- 
tween spiral arms and the surrounding material. Again, 
a similar case of competing effects comes into play, given 
that spiral arms are also invoked as a mechanism to col- 
lect inefficient, low rj^"' gas into bound, high efficiency 
objects (e.g., [^ [2000t iKoda et al.l [2001 . Similarly, 
bars may both suppress inflowing gas and concentrate 
material into nuclear starbursts fe.g.. [Sa k amoto et ahl 
[T999l:[Jogeeeran[200llSheth et al.l[200ar 
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Though competing effects cloud a simple interpreta- 
tion of the data, the path forward here remains rela- 
tively clear: systematic measurements of the internal 
conditions in the molecular gas represent a critical next 
step and our knowledge of the CO-to-H2 conversion fac- 
tor must improve to further refine our understanding of 
how T^'^ depends on physical conditions. More imme- 
diately, we need to understand which local conditions 
drive galaxy-averaged trends seen in this section. Fi- 
nally, we emphasize that while these details are essential 
to a complete understanding of star formation in galax- 
ies, they appear less critical to explain Figure [7] than the 
CO-to-H2 conversion factor and represent, in some sense, 
a factor of ~ 2 level correction to a constant rj^p . 

4.2. Local Variations in t^"' 

HERACLES resolves our targets, so that we can in- 
vestigate the dependence of the local t^I'^ on local con- 
ditions. As Figure [7] shows, galaxy averaged properties 
have a high degree of covariance. Examining local con- 
ditions may help break this degeneracy. Figure [S] plots 
the data density (blue contours) for fixed aco sls a func- 
tion of local (kpc-scale) conditions: {top left) stellar mass 
surface density, E,, {top right) H2 surface density, Smoi, 
{bottom left) dust-to-gas ratio, and {bottom right) orbital 
time. In red we show the median and standard devi- 
ation for data binned by the abscissa. These medians 
incorporate the fraction of upper and lower limits in the 
bin, which appear as percentages running along the top 
and bottom of each plot. The gray region in each plot 
indicates where the percentage of limits exceed ~ 30%, 
indicating significant issues with completeness. 

Table [7] reports rank correlations between t^'^ and lo- 
cal physical conditions. We restrict these calculations 
to the region of approximate completeness. This mini- 
mizes biases, but prevents us from probing very low sur- 
face densities. Binning and profile wor k will allow us 
to extend these analyses in future work (jSchruba et al.l 
|2011|) . Figure ini shows little systematic variation t^'^ as 
a function of and Smoi, only a weak tendency to have 
with shorter r^°' at the highest surface densities (for 

constant aco)- These trends have the opposite sense of 
those expected based on our examination of galaxy av- 
erage properties, where high mass, high I]hi+h2 systems 
showed longer depletion times than low mass systems. 

We observe a more significant relation between t^°' 
and D/G, one that becomes even stronger when limits 
are factored in. The functional form of the local trend 
in rQ°p as a function oi D/G closely matches the trend 
seen for whole galaxies. To some degree this reflects the 
fact that the D/G, unlike Ti^, or Emoi remains relatively 
constant across the disks of many of our targets (see Ap- 
pendix) . By contrast even a galaxy with very high total 
Af, or Afinoi will have a wide range of and Emoi- 

Figure [TU] shows as a function of the same local 
conditions after applying our "E = 100" conversion fac- 
tor. The correspondence between t™^ and D/G mostly 
vanishes if we adopt this conversion factor. Meanwhile 
the weak correlations of with E» and E,„oi remain 
weak. 

For fixed acOi Figure [9] does not support the idea that 



a fixed fraction of gas is converted to stars each orbital 
time. In the bottom right panel the red bins do not 
match the dashed line, which shows a fixed fraction of 
gas converted to stars per orbital time. Adopting the 
"E — 100" conversion factor changes this picture some- 
what, as the central aco depressions lead to shorter 
in regions with short orbital times. With the "E = 100" 
aco, a fixed efiiciency per orbital time becomes a rea- 
sonable description of the data below Torb ~ 100 Myr. 
Most of our data have longer Torb than this, so that 
this statement relates mostly to the inner parts of galax- 
ies. Regardless of acOi a fixed efficiency per orbital time 
does not appear to describe most of our data. Instead, 
this may be reasonable description of integrated galax- 
ies across a wide range of luminosities (jKennicuttll 1 998cl 
[Paddi et al. 2010; Gcnzcl ct al. 2010) or relevant to inner 
regions where the orbital timescale becomes comparable 
to the internal dynamical time of bound clouds. Treat- 
ment of aco makes a large difference to the results in 
this plot, highlighting the need for improved constraints 
in the inner disks of galaxies where variations in physical 
condition presumably dominate aco variations. 

In Section Hrj wc discuss the idea that the conversion of 
the diffuse, unbound ISM from atomic to molecular may 
lead to the high t"^!,"' found in massive galaxies. The local 
drivers of such trends are not immediately obvious from 
Figure[9] In Figure[TT]we directly plot r^"' as function of 
the local molecular fraction, /moi = Enioi/(Snioi + Shi). 
For /inoi > 0.5 molecular gas dominates the ISM mass 
budget. We do not know whether this molecular gas is 
organized into bound, star-forming clouds or diffuse, in- 
ert material but if high /„ioi does correspond to a higher 
fraction of diffuse molecular material, we might expect 
to observe a general increase in t^°^ as /moi increases. 
As Figure [TT] shows, our completeness severely limits this 
calculation, restricting us to /moi ^ 0-5. Above this value 
regardless of how wc treat aco we find little or no cor- 
relation of T^°' with /moi- 

Comparison of to local conditions thus reveals 
the same strong trend with D/G observed for galaxy- 
average properties but only weak trends with other pa- 
rameters, including the molecular fraction. We observe a 
suggestion of decreased t™°^ at high S* or high Emoi and, 
after apply depressed aco in galaxy centers, we find a 
weak correlation of orbital time and t™^ for short orbital 
times. 

We emphasize that these represent our broad-brush 
results. Our database will allow deeper exploration 
via detailed analysis of individual galaxies, deep Y'^'^o- 
files, varying weighting and normalization, and st ack- 
ingj Schruba et all 120 111) . Kinematic analvs is Ce.g.. iTan 
200c ) and the inclusion of outer disks (e.g., iBigiel et ah 
2OIOD should yield the lever arms to better understand 
the impact of local conditions on star formation in molec- 
ular gas. Indeed, as we discuss in the next section, sig- 
nificant peripheral evidence points to the existence of 
significant environmental dependencies of t™^. Such ef- 
fects are not immediately evident from the simple tests 
that we carry out here, however. Beyond the clear cor- 
relation of T^p' with D/G, which we interpret as likely 
due to conversion factor effects, the absence of "smoking 
gun" correlations represents the main first-order result 
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TABLE 7 

Rank Correlation of t^°' with Local Conditions 



Quantity Fixed aco ^ = 100 c«CO S = 50 oco 



S. -0.03 ±0.02 -0.01 ±0.02 -0.08 ± 0.02 

Smoi -0.06 ±0.02 -0.04 ±0.02 +0.15 ±0.02 

/moi = Smoi/{S,„oi + Shi) +0.04 ± 0.03 +0.01 ± 0.02 +0.04 ± 0.02 

Dust-to-gas ratio +0.39 + 0.02 +0.09 + 0.02 -0.23 ±0.02 

Orbital Time -0.05 ± 0.02 +0.07 + 0.02 +0.12 + 0.02 

Galactocentric Radius (r/r25) -0.05 + 0.02 +0.09 + 0.02 +0.13 + 0.02 



Note. — Rank correlation between for individual kpc-resolution lines 

of sight and local conditions for regions where we are reasonably complete. 
Uncertainties give the Icr scatter of the rank correlation about zero obtained 
when randomly repairing the data and accounting for an oversampling factor of 
4. The three columns report results for different assumed CO-to-H2 conversion 
factors (32}. We report results only for the range over which < 30% of our data 
are limits. This is approximately: S, > 40 Mq pc~^, Smoi > 4 Mq pc~^, 
/mol > 0.5, D/G > 0.006, Torb < 0.23 Gyr, and rgai < 0.55 r25. A strong 
correlation with the dust-to-gas ratio is evident for fixed aco ■ For ^ variable 
oco weak correlations. 



10.5 



38% 35% 34» 17% 10%.. ^. S% 2% 1% 0% 
Frelction of. loVer^.limlts 



10.0 



8.0 




F^Ctioti df^pper limits 
47Z3i% 19% 8% 3X i% 3% 3% 3% 3% 



1.0 1.5 2.0 2.5 3.0 

log,o Stellar Surface Density [M© pc'*] 



]fi9i 43% '88% 13%' 4% ^\0%- ■ 0% " 6% 0% 
Frecllon ot/,fov^ linnils 



8.0 




4J:ef'l6?8;.i!?..6*. 0% 0% 0% 0% 0% o% 



0.0 0.5 1.0 1.5 2.0 



2.5 




logio Dust-lo-Gas Ratio log,o Orbital Time [yr] 

Fig. 9. — Local t^°' = Smoi/EgpH for fixed oco ^ ^ function stellar mass surface density (top left), H2 surface density (top right), 

dust-to-gas ratio (bottom left), and orbital period (bottom right) for a fixed aco- Gray points show individual lines of sight and shaded 
blue contours show density of data. Red points show the median (including relevant upper and lower limits) and standard deviation in 
logiQ T^°' binned by the a;-axis. Percentages along the top and bottom indicate the fraction of lower and upper limits in each bin. Gray 
regions indicate where > 30% of our T!j^°' estimates are limits so that completeness represents a serious concern. The dashed curve in the 

lower left panel shows the expected relation for an aco-dependent D/G and the dashed line in the lower right panel shows the expectation 
for a fixed fraction of molecular converted to stars per orbital period. 
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Fig. 10. — As Figure|9]but after application of our "S = 100" conversion factor following lWolfire et alJ l|2010l . Equation[4ll. The 
factor correction removes the strongest visible trend in the data. 
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Fig. 12. — Evidence for spatial correlation among t^°' within 
galaxies. We plot RMS scatter in logj^Q t^°' (j/-axis) as a function 
of spatial resolution (x-axis) for a subset of nearby, large galaxies. 
The scatter in t^"' changes more slowly as a function of reso- 
lution than one would expect averaging uncorrelated data. Ex- 
trapolating the trend that we observe to much higher resolution 
yields rea sonable consi stency with (scarce) high resolution litera- 
ture data I IBlanc et al.l2 009: Schruba ct al. 2010). We illustrate the 
median trend, /3 = 0.5, while /3 = 1 would be expected for a disk 
filled with uncorrelated star-forming regions. This low /3 reflects 
systematic variations in T^°'on intermediate scales, either due to 

synchronization of sta r formation on intermediate scales (e.g., see 
IFeldmann et al]|20111) or systematic but still undiagnosed environ- 
ment dependencies of the efficiency of star formation. 



TABLE 8 

Scale Dependence of Scatter 



Galaxy Averaging Index, /3 



NGC 0628 


0.5 


(0.2-0.6) 


NOG 2403 


0.0 


(0.0-0.3) 


NGC 2903 


0.6 


(0.4-0.6) 


NGC 3627 


0.5 


(0.4-0.6) 


NGC 4725 


0.3 


(0.3-0.4) 


NGC 4736 


0.8 


(0.1-0.8) 


NGC 5194 


0.6 


(0.5-0.6) 


NGC 5457 


0.5 


(0.1-0.5) 


NGC 6946 


0.3 


(0.2-0.4) 



Note. — Averaging index, /3, 
for well-resolved galaxies. We 
quote the be st-fit /3, defined in 
Equation 1121 estimated from a 
minimization using Ha-|-24/^m 
and our variable "S = 100" ctco- 
In parenthesis, we give the range 
of values measured as we vary SFR 
tracer and change adopted aco- 
For comparison, we expect /3 1 
for the case of uncorrelated aver- 
aging in a thin disk. 



of this section. 



4.3. 



Spatial Correlation of t^""' 



We have examined the explicit dependence of t^'^ on 
various host galaxy properties and local physical condi- 



tions. We find many systematic dependencies on host 
galaxy, but the picture relating r^°' to local conditions 

remains more ambiguous. Considering the scatter in t"^°' 
as a function of scale offers another way to approach this 
problem. For star formation uncorrelated on the scale of 
individual regions in a two-dimensional disk, there is a 
clear analytic expectation for the functional form of this 
averaging. Deviations from this scaling can reveal the 
degree to which adjacent regions share the same t^'^, or 
at least appear synchronized. 

At small, ~ 10-50 pc, scales observations of the 
Milky Way and the nearest galaxies resolve Ssfr and 
Smoi into discrete star forming regions, clouds, and 
clusters. These individual elem ents have distinct ages 
and evolutionary sequenc es (e.g. iKawamura et al.ll2009l : 
iFukui &: Kawamrirall2010D so that the ratio of molecular 
gas to stars in a region and the emission of star for- 
mation rate tracers both evolve as a function of time 
(see discussion and plots in L12). As a result, scaling 
relations between SFR and H2, which capture the time- 
averaged relation between gas and star formation, emerge 
only after averaging together many distinc t regions (see 
iSchruba et "an[20To[ iFteldmann et al][20TTI ). 

In a disk with fixed rj^"' if these individual regions 
form stars independently then we will expect the scatter 



in r^J^p over part of the galaxy to go as vA^ , where 
A^ is the number of star-forming region in that part of 
the galaxy. For a region of extent I in a smooth disk, 
A^ oc P, so that the expectation for the "uncorrelated 
case" is <T (X l^^ . 

Deviations from this scaling will emerge if t^'^ varies 
systematically on large scales across the disk. In that 
case a high (or low) Tj^°' in one region is likely to be re- 

fiected with a similar tJ^"' in the adjacent regions. Thus, 
if Tj5^p is correlated between two adjacent regions, we ex- 
pect a weaker dependence of the scatter in r^"' on scale. 
We would expect this to occur in the case that t^"' has 
real, but still undiagnosed, dependence on local physical 
conditions. Moreover large scale dynamical effects like 
spiral density waves, bars, supernova explosions, may 
synchronize the star formation process on scales larger 
than a single cloud. 

With ~ kpc resolution, HERACLES offers limited abil- 
ity to measure the scale dependence of scatter over a 
large dynamic range, but we have identified a subset of 
large, nearby galaxies where we can measure the scat- 
ter in Tj5^°' at linear resolutions from 0.6-2.4 kpc. These 
are labeled as our "multi-scale" sample in Table [TJ We 
use this sample to measure the scatter in log^Q r^"' as a 
function of linear resolution. To do this, we convolve each 
galaxy in our "multi-scale" sample (Table [1]) to have lin- 
ear resolution 0.6-2.4 kpc and measure the RMS scatter 
in logj^o '''Sep across the galaxy at each resolution. Doing 
so, we make no correction for inclination, so that this 
exercise consists of placing targets at larger and larger 
distances. 

Figure [12] plots the results of this exercise for each re- 
solved galaxy. We show cr, the RMS scatter in log^p T^lp 
as a function of spatial resolution. The figure shows re- 
sults calculated using the "S — 100" acoj which removes 
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a significant internal gradient in r^°' from NGC 5457. 
For comparison we plot the scatter in r^"' measured at 
high resol ution in M3g| bylSchruba et"all (|20ia stars) 
and M51 (iBlanc et al.ll2009l diamonds). 

Figure [T2] shows a steady increase in scatter with im- 
proving linear resolution. We characterize this scale de- 
pendence of the scatter via a power law, 



C (0 = '''600 



I 



600 pc 



(12) 



where I represents the spatial resolution, crgoo is the scat- 
ter in T^l'^ at 600 pc resolution and the power-law index 
/3 measures the rate at which changing the resolution 
changes the measured scatter in r^"' . We report the 
best-fit /? for each multiscale target in Table [51 As de- 
scribed above, we expect (3 — 1 for uncorrelated star 
formation in a disk. 

In most cases the best-fit averaging index, f3, is ~ 0.5, 
significantly less than the f3 — 1 expected for uncorre- 
lated, fixed-efficiency star formation. Based on the pre- 
vious sections, we expect that this refiect real r^""' varia- 
tions or correlated systematic uncertainties in our phys- 
ical parameter estimation (e.g., undiagnosed aco varia- 
tions). That is, this is another way to see the subtle but 
real systematic variations in t^°' considered in the last 
two sections. Alternatively, Figure [T^ and Table [5] could 
reflect a high degree of synchronization, with adjacent re- 
gions likely to be at the same stage of the star formation 
process and thus show similar ratios of star formation 
tracers to CO emission. Either synchronization or real 
r^p' variations might be achieved by dynamical phenom- 
ena at small scales and systematic efficiency variations 
may arise from dependence on local conditions that we 
have yet to identify. 

This weak, f3 ~ 0.5, scale dependence agrees 
with predictions base d on numerical simulations by 
iFeldmann et al.l ()201lD , who discuss /3 in terms of dimen- 
sionality. As an example, /? = 0.5 would be consistent 
with the degree of correlation expected by spiral arms 
or any other phenomenon that synchronizes star forma- 
tion along one dimens ion, though in the simulations of 
IFeldmann et al.l ()201lD it arises more generally. The low 
/3 also agrees qualitatively with numerous observations 
of highly structured star formation on scales < 1 kpc in 
nearby galaxies" for example, see the recent synthesis by 
IFlmegreen {201%. 

From the perspective of galactic-scale star formation, 
the key point from this calculation is that the scatter 
in at kpc scales depends on processes operating at 
scales larger than that of individual star-forming regions. 
That is, key information on the distribution of star for- 
mation in galaxies remains to be extracted from compar- 
ison of maps of SsFR and Smoi- This represents another 
manifestation of the overall theme of this section, that 
real second order variations in t^"' do appear visible 
in our sample. Future investigation of HERACLES and 
similar surveys will allow tests of the degree to which this 
correlation can be attributed to systematic variations in 

We infer the scatter from their measurements of CO-to- 
Ha-|-24/im near Ho peaks and CO peaks. 



■^dcp ^ function of yet-unexplored local conditions or 
to physical parameter estimation. Simultaneously, obser- 
vations with high spatial dynamic range will allow more 
a detailed diagnosis of the scale-dependence of r"^"'. 

4.4. Enhanced Efficiency in Galaxy Centers 

Our sample does not contain any true galaxy-wide star- 
bursts, but many of our targets do host n uclear c oncen- 
trations of star formation and gas (HclfereLaUSOOS), 
so that the inner kpc of our targets represents a sub- 
sample intermediate between normal disk galaxies and 
starbursts. These regions represent only a small fraction 
of the area in our targets, so exert a negligible impact 
on the ensemble of data seen in Figure [1] but they probe 
an important part of parameter space — high Smoi, high 
^SFR — and so we explicitly consider them in Figure [T51 
We plot the enhancement in r^"', calculated by dividing 

T^"p' for each point by the average for that galaxy, and 
mark points from the inner kpc of our targets with large 
symbols: blue dots for galaxies with sonie spe ctroscopic 
indication of an AGN (jMoustakas et al.ll2010l or NED) 
and green stars show star-formation dominated regions. 

The top left panel of Figure [T3] shows enhancement in 
'''dop' for a fixed aco ■ We find systematically lower t^'^ in 
the centers of our targets, both AGN ai id starbursts. We 
find this shorter t^'^ using a fixed aco ■ ISandstrom et al.l 
(j2012l ) find evidence for systematically lower conversion 
factors in the central parts of our sample. The top right 
panel of Figure [T3] shows the results of applying these 
central corrections, as part of our "S — 100" conversion 
factor. T^p' become even shorter in the central regions 
compared to the disks. Figure [T3] thus exhibits one of 
the clearest systematic effects in our sample: the central 
regions of our targets tend to have significantly shorter 
r^p' than the ensemble of disk regions. Central r^^"' are 

0.8 times the disk value with ±0.2 dex {la scatter) for 
fixed aco- The median central-to-disk ratio drops to 0.6 
with ±0.35 dex scatter our "E — 100"' case. 

This apparently real shortening of r^"' coincides with 
an increase in the C0(2-l) to CO(l-O) ratio, indicative 
of more excited gas. The bottom right panel shows that 
this lower t™°' coincides with higher C0(2-l) to CO(l-O) 

ratios (calculated with comparison to iKuno et"aLll2007D . 
The central regions of our targets show systematically 
enhanced CO(2-1)/CO(1-0) compared to the disks. This 
excitement presumably reflects the same changing phys- 
ical cond itions that drive the lower conversion factors 
found bv ISandstrom et all (I1Q12:;), underscoring that nu- 
clear gas concentrations represent a distinct physical 
regime from galaxy disks. Molecular gas in these regions 
gives off more CO emission, appears more excited, and 
forms stars more rapidly than molecular gas further out 
in the disks of galaxies. 

Not all nuclear regions exhibit shorter t^I'^. A suffi- 
cient but not necessary condition to find such enhance- 
ments in our data is: (1) to lie within the central kpc of 
a target and (2) have C0(2-l) intensity > 15 K kms~^. 
We indicate this criteria using a red line in the bottom 
left panel of Figure [T31 This does not precisely equate 
to a Emoi threshold because of ambiguities introduced by 
aco; the figure shows that central parts of galaxies with 
bright CO emission tend to show lower . 
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Fig. 13. — Enhanced efficiency in galaxy centers. We plot t^°' normalized to the galaxy average as a function of galactocentric radius 
for a fixed oco (top left) and a "S = 100" conversion factor (top right). We also plot enhancement in t^°' as a function of CO (2-1) 
intensity (bottom left) and show the enhancement in the CO(2-l) to CO(l-O) line ratio relative to the disk average as a function of radius. 
Individual gray points show kpc-resolution lines of sight; contours show data density; black-and-white points and error bars show median 
and RMS scatter in binned data. Blue points and green stars show data from the inner kpc of our targets — blue points show systems with 
a spectroscopic classification indicating the likely presence of an AGN tMoustakas et al..,2010. ) , green points show star formation dominated 
nuclei. Galaxy centers show shorter depletion times and enhanced line ratios relative to galaxy disks, indicative of more excited, more 
efficiently star-forming gas. The effect appears even stronger once we account for variations in oco (Sandstorm et al.) and essentially all 
CO- bright nuclear regions show some level of enhancement. 



IDaddi et all (|2010l) and IGenzel et all (j2010t ) argue for 
multiple "modes" or "sequences" of star formation, in 
which the shorter dynamical time and higher density in 
starburst galaxies lead to shorter depletion times at fixed 
Smoi- Figure [T3l supports this idea inasmuch as it shows 
that physical conditions other than gas surface density 
play an important role determining t^'^. The popula- 
tion of low T^p' in Figure [12] appear to be set by factors 
other than kpc-scale Smoi alone: first, identifying these 
lines of sight requires knowing that the points line in the 
centers of our targets; second, t^"' exhibits a wide range, 
presumably set by other physical parameters. Figure [13] 
does not support the idea of two cleanly distinguished 



sequences, suggesting instead that a continuum of 
exist in the central regions of our targets; we plot the 
histogram of central t^'^ in Section [5] 

What drives these enhancements? Beam dilution cer- 
tainly plays some role. Our kpc resolution will aver- 
age out nuclear gas concentrations with surface densities 
well in excess o f the disks of our targets but small spa- 
tial extent (e.g.. lJogee et al.l[2005[ ). We expect that high 
pressure, driven by the deep potential well in the cen- 
tral parts of galaxies, also plays a role driving gas to 
higher densities. This effect is seen in our own Galaxy 
dOk a et al 2001} and others (Rosolo wskv & Blitz 2005). 
iDaddi et all (|2010h and lGenzel et aO (|2010t ) suggest that 
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the shorter dynamical timescales in starburst galaxies 
may also play a key role and our data may offer some 
tentative support for this (Figure[TO]). For this paper, the 
key result is that the shorter t^I*^ do exist, so that the 
inner parts of disk galaxies represent a kind of "transi- 
tion regime" between the disks of spirals and galaxy-wide 
starbursts. Follow-up interferometry and spectroscopy 
contrasting physical conditions in the nuclear regions 
with those in disks will yield more insight. 

5. DISCUSSION AND CONCLUSIONS 

We study the relationship between molecular gas sur- 
face density, Emoi, as traced by CO (2-1) emission and 
SFR surface density, Ssfr — traced by combinations of 
Ha, FUV, and IR emission — at 1 kpc resolution across 
the optical disks of 30 nearby spiral galaxies. Broadly, we 
demonstrate two conclusions, which we illustrate using 
histograms of t^'^ in Figure [TH a first-order simple cor- 
respondence between Smoi and Ssfr and second-order 
systematic variations in the apparent molecular gas de- 
pletion time, Tj^°p', including lower values in nuclear star- 
bursts and low, mass low-metallicity galaxies. Some, but 
not all, of these variations may be explained by invoking 
a dependence of the CO-to-H2 conversion factor, aco, 
on dust-to-gas ratio 

First, molecular gas and star formation are tightly cor- 
related in both individual galaxies and over individual 
kpc-resolution lines of sight. Their ratio, the molecular 
gas depletion time, r^"' appears to first-order constant 
across the disks of massive, large star-forming galaxies, 
with a median value t^'^ = 2.2 Gyr, a scatter of w 0.3 
dex, and systematic uncertainty ~ 60% (see the nearly 
lognormal histogram in the left panels of Figure [T4| . We 
arrive at this conclusion testing a wide range of method- 
ologies to trace the recent star formation rate and distri- 
bution of molecular gas. This includes the interchange 
of Ha and FUV emission, varying treatments of the in- 
frared "cirrus", substitution of literature CO (1-0) for 
our CO (2-1) maps, and adoption of variable CO-to-H2 
conversion factors. We demonstrate that our measure- 
ments agree with a large collection of literature data 
from the last decade fSection 13.41) . with all data occu- 
pying a common region of Smoi-SspR parameter space. 
Our data reinforce and extend a consensus for a "large 
disk galaxy" value of r^"' « 2 Gyr for matched assump- 
tions about the CO-to-H2 conversion factor and stellar 
IMF with « 0.3 dex scatter. 

Adopting a forward -modeling approach similar to that 
of lBlanc et all (pOOl . we derive a best-fit power law in- 
dex of iV « 1 ± 0.2 for EsFR oc E^^j (Section |33]) for 
a fixed aco and show that individual galaxies exhibit 
a distribution of N with N mostly in the range 0.8-1.2. 
However, we stress the inadequacy of a power law to cap- 
ture important changes in physical conditions other than 
Smoi and the substantial uncertainty in suc h fits. Rein- 
forcing the conclusions of lBlanc et al.l ()2009[ ). we caution 
that the commonly used combination of "sigma-clipping" 
and bivariate fitting has the potential to substantially 
bias results (see Appendix), yielding seemingly discor- 
dant values of N even when the data substantially agree. 

Our second major conclusion is that with a broad sam- 
ple spanning a wide range of physical conditions, system- 
atic variations in the apparent t^I°^ emerge both among 



and within galaxies (see the width of the distribution in 
the top middle panel of Figure [14]) . We show systematic 
variations of galaxy-average apparent t^"^ as a function 
of many host galaxy properties: stellar mass, rotation 
velocity, metallicity, dust-to-gas ratio, average gas sur- 
face density, and morphology. These variations have the 
sense that low mass, low metallicity, late type galaxies ex- 
hibit shorter apparent r"^°' than high mass galaxies. The 
trends persist, though weaker, even into the high mass 
region, A/» > 10^° M©, and agree well with those seen 
in the COLDGASS sample (jSaintonge et al.l[20Tl . We 
emphasize "apparent" because these variations appear 
to be a mixture of real changes in the rate at which gas 
forms stars and biases in physical parameter estimation. 
Adopting a CO-to-H2 conversion factor that depends on 
the dust-to-gas ratio can explain many of the strongest 
variations with host galaxy properties; note the narrow- 
ing from the top middle panel of Figure fT4l to the bottom 
middle panel. Our ability to examine residual trends re- 
mains restricted by the limited precision with which we 
know aco J but our best estimate is that correlations do 
remain between the real (aco-adjusted) t^'^ and galaxy 
mass, average gas surface density, and perhaps several 
other quantities. 

We also examine how t^I'^ varies as a function of local 
conditions. We find two strong relationships: the appar- 
ent r^p' calculated for a fixed aco varies systematically 
as a function of dust-to-gas ratio and we observe system- 
atically lower and widely varying t^I'^ in the inner kpc 
of our targets (see the top right panel of Figure [M] We 
interpret the first as indicating important variations in 
aco and, as with the galaxy-integrated case, show that 
application of a dust-to-gas ratio dependent aco can ex- 
plain much of the observed trend. The lower t^I'^ in 
galaxy centers appears real and robust to aco consider- 
ations. Indeed, iSandstrom et al.l (|2012| ) find low aco in 
the central parts of many of our targets, which implies 
even lower t^"'. The resulting t^^'^ in the inner parts 
of galaxies varies widely (see the right panels in Figure 
fT4|) . providing strong evidence that environmental fac- 
tors do drive rj^"', and thus the star formation rate, in 
these regions. 

Strong local trends beyond those linking t^I'^ to the 
dust-to-gas ratio and nuclear starbursts elude our present 
analysis, though the correlation of t^"' with integrated 
properties and weak dependence of scatter on scale 
strongly suggest their presence. We suspect that the 
limited resolution and remaining imprecision in physi- 
cal parameter estimation so far obscure these trends in 
HERACLES. Comparison to existing Herschel data and 
spectroscopy of other molecular lines should improve our 
ability to understand the origins of the physical scatter 
in Tj^^p'. Perhaps just as important, at kpc resolution 
significant averaging has already occurred, especially in 
regions of high Ssfr- The star- forming ISM hosts many 
competing effects: for example, shear may both suppress 
collapse and lead to more frequent cloud collisions; high 
pressures may lead to both denser clouds and a substan- 
tial diffuse molecular ISM; spiral arms may both collect 
material and suppress collapse via streaming motions. 
At our kpc resolution, such trends will be subtle as com- 
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Fig. 14. — Distribution of molecular gas depletion time, t^'^, in HERACLES. We show histograms for (left) all lines of sight at kpc 
resolution, (middle) whole galaxies, and {right) lines of sight in the central kpc of our targets. The top row shows results for fixed aco 
while the bottom row applies our "S = 100" conversion factor, which accounts for variable dust shielding and central aco depressions. The 
distribution of t^'^ treating lines of sight equally appears lognormal, with a depletion time 1-3 Gyr, depending on the adopted conversion 
factor, and < 0.3 dex scatter, though low-Tg'°p excursions from the lognormal distribution are already visible in this panel. The distribution 
of T^°p for whole galaxies shows more scatter for a fixed aco > reflecting systematic trends in which low-mass, low-metallicity, high sSFR 
galaxies show faint CO relative to SFR tracers. After accounting for varying ocOi there is less scatter among galaxies but some spread 
remains. Meanwhile the inner kpc of our targets shows a much broader, much less normal distribution, with both low and high t^°' 
common. Correcting for the central «co depressions observed bv Sandstrom et al.l 120121) and expected from observed line ratio variations 
exacerbates this effect, but it is clearly present even adopting fixed aco- Thus, t^°' appears constant to first order but clear systematic 
second-order variations do emerge as a function of local and galaxy scale conditions in a broad sample. 



peting effects occur inside a single resolution element. In 
the longer term, higher resolution observations of a di- 
verse sample of galaxies will be needed to diagnose the 
impact of dynamical effects on star formation. 

These results fit into a broad picture of star formation 
in galaxies as follows. Within the disks of nearby galax- 
ies, we find recent star formation correlat ed with molec- 
ular, rather than atomic or total, gas (jSchruba et al.l 
120 lit ) . The ratio of CO emission to recent star formation 
appears roughly constant within massive , star-forming 
disk galaxies (Section [51 Bi giel et al.|[20TTI ) but when ex- 
amined more closely, significant variations do emerge be- 
tween this ratio and galaxy mass ( Saintongc ct al. 2012, 
Section |4]) and dust-to-gas ratio or metallicity (Section 
m. The interpretation of these trends depends critically 
on the behavior of the CO-to-H2 conversion factor. We 
show that current best estimates may explain many of 
the observed trends in the CO-to-SFR ratio (see also 
iBolatto et all 120111: ISchruba et al.l 120121. However, we 
caution that significant work is still needed to bring con- 
version factor estimates to the precision needed to confi- 
dently interpret these trends. Our best estimate is that 
correlations in which the true molecular gas depletion 
time, T^p', increases with increasing galaxy mass or gas 



surface density do persist after accounting for aco ef- 
fects, but that these are comparatively weak. A sensible 
explanation for such trends is the emergence of a diffuse, 
unbound molecular medium at high gas surface densi- 
ties and high pressures, but direct evidence relating a 
local high molecular fraction to lower t^'^ is weak in our 
present data set. We do find good evidence for system- 
atically low r^p' in galaxy centers, with a wide variation 

in the factor by which t^'^ falls below the average for 
that disk. This supports the idea that in environments 
with high surface densities and short dynamical or or- 
bital times, environmental factors may drive t^'^ to a 
wide range of values at fixed average gas surface den- 
sity. This is a more general formulation than the idea 
of distinct "disk" and "starburst" s equences, but quali - 
tati vely agrees wi t h the picture from lDaddi et al.l ()2010f ) 
and lGenzel et al.l (|2010D . with the nuclei of disk galaxies 
occupying a regime intermediate between quiescent disks 
and merger-induced starbursts. 
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Fig. 15. — Line ratios and dust-to-gas rati os in HERA CLES, (left ) CO (2-1) / CO (1-0) lin e ratio comparing HERACLES integrated 
measurements to literature CO (1-0) surveys: IKuno et al.l (2007), Hei fer et al.l II2003I ) . [Young eFa l. (1995) and measurements for individual 
pointings obtained for foUowr-up spectroscopy by Usero et al. (in prep.). The ensemble of measurements has median 0.67 with scatter 
0.16 dex. (right) Dust-to-gas ratio profiles derived comparing IR SED modeling, HERACLES, and H I data. We plot the dust-to-gas ratios 
derived fixed conversion factor in radial profile. 



APPENDIX 
MEAN LINE RATIO FOR HERACLES 

HERACLES surveyed the C0(2-l) transition. Our adopted "fixed" aco = 4.35 has been derived considering mainly 
CO(l-O) observations and most of the literature comparing recent star formation and molecular gas considers CO(l-O). 
Rosolowsky et al. (in prep.) present a thorough analysis of the line ratio as a function of local conditions in HERACLES 
and explore the implications for aco and ■ In this study we adopt a fiducial CO(2-1)/CO(1-0) ratio of 0.7. The 
left panel in Figure lAl shows a histogram of CO(2-1)/CO(1-0) values measured for the HERAC LES survey. W e plot 
ratios derived from com paring integrated H ERACLES fluxes to those from the CO(l-O) surveys bv lYoung et all ([1995), 
iHelfer et al.l ()2003[ ). and lKuno et alj ()2007f ). We also show the results of pointed single-pixel spectroscopy carried out 
with the IRAM 30-m, part of a large spectroscopic database present by Usero et al. (in preparation). These have been 
reduced using CLASS in a standard way and aperture-corrected to match beam areas using the C0(2-l) distribution 
in HERACLES. The ensemble of measurements has median 0.67, with the ratio taken in brightness temperature units, 
and a scatter of 0.16 dex, sa 40%. Some of this scatter will be due to calibration uncertainties in the CO(l-O) data. 
However, even comparing each data set to HERACLES separately, we find internal scatter in CO(2-1)/CO(1-0) of 
20-40%. We discuss the most obvious environmental dependence of this ratio, the enhancement in galaxy centers in 
Section (see also lLerov et al.ll2009[ ). Based on Figure |X1 we adopt 0.7 as a typical line ratio. In this paper, the 
primary application of this value is to apply a "standard" Milky Way CO(l-O) conversion factor to the HERACLES 
C0(2-l) data. 

DUST-TO-GAS RATIO AND CONVERSION FACTOR CALCULATIONS 

The right panel in Figure [Al plots the dust-to-gas ratios, D/G, that we derive from a fixed aco conversion factor 
and our Spitzer SED modeling. 

Correction Due to CO-Dark Gas 

We use the recent theoretical work by iWolfire et al.l (|2010[ ) to estimate the fraction of "CO-dark" molecular gas. 
This gas lies in in regions where car bon is mostly associat ed with C II rather than CO, and so will no t be readily 
traced by maps of CO emission (see iKrumholz et al.ll20lil for a similar approach). IWolfire et al.l (|2010D present an 
expression for the fraction of mass in this "CO-dark" phase, 

/co-da.k = l-exp(^-_^j (Bl) 

where AAv.dg is the depth that the CO-dark phase extends into the cloud, measured in units of visual extinction, 
and Ay is mean extinction through the whole cloud. Ay simply depends on the product of cloud surface density — 
or equivalently mean extinction through the cloud at solar metallicity. Ay — and the dust-to-gas ratio. Ay {D/G') = 
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Fig. 16. — Distribution of best-fit parameters for 1, 000 Monte Carlo iterations designed to test the impact of observation uncertainties, 
sample definition, and fitting technique on the best-fit power law. We plot the distributions of slope n, coefficient a, and intrinsic scatter s 
across all iterations for our best SsFR estimate (H«+24/xm) and a fixed Qco-- 



A^y D/G'. IWolfire et all (poToH give an expression for lS.Ay £,Q that depends weakly on the density and metaUicity of 
the cloud 
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Where Gq is the radiation field relative to the Solar Neighborhood value, ric is the density of the cloud, and Z' is 
the metallicity (or dust-to-gas ratio) relative to the Galactic value. We neglect the second term, setting Gq = 1 and 
lie = 1 cm~'^ and use the D/G' in place of Z' . /S.Av^dg does not vary much with any of these quantities, so that 
/co-dark mainly depends on the extinction through the cloud, Av{D / G'). 
We scale the conversion factor by this expression, so that: 
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Adopting the conversion between dust and column density of iBohlin et al.l ()1978[ ) and a fiducial average molecular 
cloud surface density of 100 pc^^, so that Ay ~ 4.9, we calculate the correction for "CO-dark" gas via: 



cco-dark {D/G')^aco {D/G') /aco {D/G' = 1) 

= aco {D/G') /aco {D/G'^1) 

( 4.0 X 0.53 \ , /4.0 X 0.53 
= exp — — — — / exp 
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where Eioo is the assumed universal surface density of molecular clouds normalized to 100 pc ^ and D/G' is the 
dust-to- gas ratio normalized to the Milky Way value. This is Equation [5l 

Calculation of Dust-to- Gas Ratio and aco From Profiles 
ISSUES IN POWER LAW FITTING 

Monte Carlo Uncertainty Estimates 

If we treat our goodness-of-fit statistic (Equation [11]) as then the implied statistical uncertainties on our fits are 
very small (the + 1 surface implies uncertainties order 1%, as in iBlanc et al.ll2009D . At some level, this is accurate: 
given our fitting approach and data set, the best-fit power law is heavily constrained. However, this does not refiect 
our real best estimate of the uncertainty in the underlying relationship between Ssfr and Smoi- To estimate a more 
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TABLE 9 

Monte Carlo Estimates of Uncertainty Power, Law Fit Parameters 



Source of Uncertainty 


Scatter in n 


Scatter in logjQ a 


Scatter in s 


Statistical Noise 


0.05 


0.01 


0.02 


Calibration Uncertainties 


0.12 


0.06 


0.04 


Robustness to Removal of Individual Data 


0.07 


0.03 


0.02 


Robustness to Removal of Galaxies 


0.10 


0.06 


0.05 


Choice of Fitting Parameters 


0.06 


0.02 


0.02 


Overall Uncertainty 


0.15 


0.09 


0.05 



Note. — Each entry gives the Itr standard deviation in the parameter across 100 Monte 
Carlo iterations (1,000 for the overall uncertainty). We quote results for our best-estimate 
^SFRi Ha-|-24/im corrected for IR cirrus, and a fixed aco- 



realistic uncertainty, we carry out a series of Monte Carlo simulations. We begin with our measurements, estimates of 
Smoi and SsFR at each point. We then examine the effects of: 

• Statistical noise. For each point, we add normally-distributed noise of the appropriate magnitude to each Emoi 
and log-normally distributed noise of magnitude 0.15 dex to Ssfr (see Section^. 

• Calibration uncertainties. Calibration issues include uncertainties in the overall flux scale of the data or uncer- 
tainties in the conversion to physical parameters, e.g., due to variations in CO line ratios, dust properties, or 
stellar populations. In each case, these wi ll tend to op erate galaxy-by-galaxy. We take these to be log-normally 
distributed with a magnitude of 0.15 dex (|Lerov et al.ii200i^ L12). We scale all of the data for each axis in each 
galaxy by a single, randomly generated factor. 

• Robustness to removal of individual data. We test the robustness of the results to the removal of individual 
data using a standard bootstrapping approach. We resample the data, allowing repeats, to produce a data set 
matched in size to the original. 

• Robustness to removal of galaxies. We also test the robustness of our fit to the removal of a whole galaxies using 
a bootstrapping approach. Instead of resampling the ensemble of measurements, we resample the list of galaxies, 
allowing repeats, to produce a sample of 30 (non-unique) galaxies. 

• Choice of fitting parameters. A subtle point in our fitting is how to handle data in cells with low (or zero) model 
probability. A sin gle outlying daturn can dramatically skew the results if this issue is not properly treated. By 
default, we follow iBlanc et all (|2009L and G. Blanc priv. comm.) and never allowing the weight for a single cell 
to be lower than expected for N^^^^-^ = 1 (i.e., we cap the denominator in Equation [TT] at this level). We test the 
impact of this choice and our adopted grid cell size (0.125 dex). We allow the bin size and the minimum weight 
per cell to vary by up to a factor of two either higher or lower, with equal probability across the range. 

We test each of these sources alone and report the results for Ssfr traced by Ha and 24/xm and Smoi calculated from a 
fixed aco in Table IHl Figure [T51 shows the distribution of fitted parameters across all runs for our best Ssfr estimate. 
The strongest contributors to the overall uncertainty are galaxy-to-galaxy calibration uncertainties or sample definition. 
Our large data set and good S/N makes statistical noise over individual lines or sight or removal of individual data 
minor concerns. Choice of fitting parameters do impact the overall results, but do not dominate the uncertainty. 

In the end, we derive our overall uncertainty on the fits from the scatter in the best fit parameters across 1,000 
Monte Carlo iterations that include all of these effects. Because we can only add noise or decrease our sample size, 
each Monte Carlo iteration operates on a data set of inferior quality to that used for our best estimates (Table S]). We 
consider these uncertainties realistic but conservative. 

Power-Law Index Biases in Bivariate Fits to Clipped Data 

Fitting a relationship betwee n Ssfr and Smoi entails conducting a bivariate fit to data with limited sensitivity. The 
approach of IBlanc et al.l ()2009( ) incorporates upper limits. However, it has been common practice to clip data at some 
signal-to-noise threshold in CO (Smoi) and then t o carry out a bivariate linear fit relating logi n Ss fr to logj^Q Smoi 
above this threshold, e.g., using the OLS bisector (jlsobe et al.lll990[ ) or FITEXY ([Press et al.lll992h methods. This 
clip-and-fit approach can introduce a significant bias into the fit power law index. We note this effect and its magnitude 
here. 

The bias arises from the interaction of the clipping with fitting techniques that attempt to minimize the two- 
dimensional distance between a point and the fit line (in log- log space) . Consider data that are intrinsically defined by 
a power law, a line in log-log space, but exhibit significant scatter orthogonal to the line. Now consider data near the 
clipping threshold in Emoi, the quantity that typically represents the limiting observable. Data that scatter to high 
EsFR and low Emoi will be discarded from the analysis due to the clipping. Data that scatter to low Esfr and high 
Smoi will be included in the analysis. Preferentially adding data at low Ssfr and high Emoi will tend to skew the fit 
towards steeper slopes. The left panel in Figure [TTl illustrates the effect for simulated data. 
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Fig. 17. — (left) Example of the biases introduced by clipping data on one axis (e.g., due to limited S/N) before carrying out a bivariate 
fit. Points show a simulated data set in which SgpR oc Smol by construction (red line) with 0.3 dex scatter introduced to each axis. We 
clip data below Sniol = 20 Mq pc~^ so that gray points are discarded and black points are fit. We then fit the black points using the OLS 
bisector method and plot the result in blue. The result is a bias towards higher slope, (right) Derived power law index from using the 
OLS bisector to fit noisy, low dynamicrange data with clipping applied to one axis. The underlying model has an intrinsic n = 1.0 slope 
with Emol drawn from our data and log-normal scatter of the indicated magnitude applied to each axis. The grayscale and contours report 
the best-fit power law index after the data are clipped at the x-axis value in one axis and fit using the OLS bisector technique. Biases of 
several tenths up to as much as four tenths are possible within the range of values adopted by recent studies. 

We estimate the magnitude of the effect using a Monte Carlo simulation. We consider samples of 10,000 data points 
with Emoi randomly drawn from our data set. We assume an underlying linear relation, Ssfr oc Smoi, so that the 
true n ~ I. We set r^""' — 1 Gyr, but the choice is arbitrary. We realize 100 samples of 10,000 data points each. 
For each sample we introduce log-normal scatter of equal magnitude to each axis and then clip the data at a series of 
5^moi values from to 40 M0 pc^^. We then carry out an OLS bisector fit to the data and record the best-fit power 
law index. The right panel in Figure [T7] plots the average fit index across all Monte Carlo iterations as a function of 
the scatter introduced into the data (y-axis) and the threshold imposed (x-axis). We find that for commonly used 
thresholds (Emoi 10-20 M© pc~^) and typical observed scatters (~ 0.3 dex) a significant bias can be introduced to 
the best fit index, often shifting it from its true value of 1.0 (by construction) to ~ 1.3 or 1.4. 

Note that Figure [17] illustrates the problem but that the exact magnitude of the bias will depend on the dynamic 
range, noise, and underlying relationship in the data set studied. Also note that conversely, clipping based on Ssfr 
will tend to bias the index in the opposite direction, to lower n. Clipping based on both Ssfr and Smoi may lead 
to offsetting biases but does not represent a proper substitute for a rigorous treatment. Finally, we emphasize that 
although we do clip data below Emoi = 5 Mq pc~^ in our fitting (Section I3.3P the bias describe here will not affect 
our results because we model the data distribution rather than carry out a bivariate fit. 

SsPR - Emol relations for individual galaxies 

Combinin g ou r measurements into a single data set obscures real differences among galaxies. Table [2 Figure [3 
and Section 14.11 show systematic variations in average r^"' and power-law index among galaxies. In Figures [18] and 
[TO] we plot SsFR as a function of Emoi for individual galaxies. In the background, we plot contours for the combined 
distribution shown in Figure [T] We show results for our best SFR tracer, Ha-f24/im, and a fixed CO-to-Il2 conversion 
factor. Points from the central part of the galaxy, Tgai < 0.1 appear in red. Points from the rest of the galaxy 
appear in gray. 
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Fig. 18. — SgFR, estimated from Ha+24/im, as a function of Snioli estimated from HERACLES CO(2-l) data and a fixed «cOi for 
individual galaxies. Data from the designated galaxy appear in gray, with data from the inner 0.1 r25 marked in red. In the background, 
we plot data density contours from the ensemble of all measurements (Figure [l} . Annotation and symbols are otherwise as Figure [l] 
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